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Abstract We explore the detailed polarization behaviour of pulsar 
0823 + 26 using the technique of constructing partial ‘mode-separated’ 
profiles corresponding to the primary and secondary polarization modes 
The charactenstics of the two polarization modes m this pulsar are 
particularly interesting, both because they are anything but orthogonal and 
because the secondary mode exhibits a structure seen neither in the primary 
mode nor in the total profile The new leading and trailing features m the 
secondary mode, which appear to represent a conal component pair, are 
interpreted geometrically on the basis of their width and the associated 
polarization-angle traverse as an outer cone 
If the secondary-mode features are, indeed, an outer cone, then ques- 
tions about the significance of the pulsar’s postcursor component become 
more pressing. It seems that 0823 -f- 26 has a very nearly equatonal 
geometry, m that both magnetic poles and the sightline all fall close to 
the rotational equator of the star We thus associate the postcursor 
component with emission along those bundles of field lines which are also 
equatorial and which continue to have a tangent in the direction of our 
sight hne for a sigmficant portion of the star’s rotation cycle. It seems 
that in all pulsars with postcursor components, this emission follows the 
core component, and all may thus have equatorial emission geometnes No 
pulsars with precursors in this sense - mcludmg the Crab pulsar -are 
known. 

The distribution of power between the primary and secondary modes is very 
similar at both 430 and 1400 MHz. Our analysis shows that in this pulsar 
considerable depolarization must be occurrmg on time scales that are short 
compared to the time resolution of our observations, which is here some 
05-10 milliseconds. One of the most interesting features of the mode- 
separated partial profiles is a phase offset between the primary and secondary 
modes. The secondary-mode ‘mam pulse’ arrives some 1.5 ± 0.1° before the 
primary-mode one at 430 MHz and some L3 + 0.1° at 21 cm Given that the 
polar cap has an angular diameter of 3.36°, we consider whether this is 
a geometric effect or an effect of differential propagation of the two modes in 
the inner magnetosphere of the pulsar 

Key words. Pulsars — polanzation — PSR B0823 -I- 26. 
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1. IntroductioD 

Pulsar B0823 + 26 is one of the oldest known pulsars, being the first or second pulsar 
discovered at the Arecibo Observatory (Craft, Lovelace & Sutton 1968) With a penod of 
531ms and a spmdown of 172 x 10"^^ sec/sec -yielding a value of 3.44 -its 

timing properties could not be more ordinary. However, Backer et al (1973) discovered 
that the pulsar has both an interpulse (IP) and a postcursor (PC) component The pulsar is 
also relatively close (DM 19.5pc/cm^), bnght, and exhibits a turnover in its radio- 
frequency spectrum only below about 50 MHz (Sieber 1973, Izvekova et al 1981) For all of 
these reasons, it has been possible to study it over a remarkably wide mterval of frequency 
On the high end it has teen detected at Bonn by Bartel et al (1978) at 14 8 GHz, and 
observations on the low end have teen reported at 26 MHz (Philhps & Wolszczan 1992) 
Over this enormous frequency range the form of the pulsar’s profile changes 
remarkably little [see Hankins & Rickett (1986) for a senes of profiles spanning 135 
to 2380 MHz] The MP exhibits a nearly Gaussian-shaped profile at both the very highest 
and lowest frequencies, only distorted by what is probably interstellar scattenng at 
26 MHz Its width does change, however, from temg practically constant above 1 GHz, to 
growing with an asymptotic dependence of about /"° below about 300 MHz (Ranbn 
1983b, hereafter Paper II) The half-power width of the MP, mterpolated to IGHz, was 
determined in Rankin (1990, hereafter Paper IV) to be 3.38°, very close to the 2457P^^^ 
angular diameter of the pulsar’s polar cap at the surface of the star 
The MP-IP spacmg is very close to 180° and could be independent of frequency, 
although a slight decrease with frequency is also compatible with the observations 
(Hankins & Fowler 1986). The MP-PC separation, however, increases monotomcally with 
frequency between about 400 and 2400 MHz as (Hankms 8c Fowler 1986), and 
Arecibo 130-MHz observations further confirm this behaviour (Hankins & Rankin 1994) 
The relative intensity of the IP is strongest at 400 MHz, with a peak intensity of just under 
1 per cent that of the MP, and decreases steadily at both higher and lower frequencies The 
PC, in turn, is somewhat stronger than the IP, being about 3 per cent of the MP at 
400 MHz At frequenaes above about 400 MHz the relative intensities of the PC and IP fall 
off at about the same rate ( /*" ” M Below 400 MHz, the relative intensity of the IP decreases 
very rapidly, whereas that of the PC mcreases, at least to about 100 MHz, though far less 
steeply than at higher frequencies (Hankins & Rankin 1994). 

Full penod, average polanmetry has been published for pulsar 0823 + 26 at 430 MHz 
(Rankin & Benson 1981), 1400 MHz (Rankm etal 1989, Blaskiewicz etal 1991), and 
10 55 GHz (Xilouns et al 1995), all but the last of which also mcludes some information 
about the PC and/or the IP Individual-pulse polarimetnc studies have been earned out at 
430 MHz (Rankm et al 1974; Backer & Rankm 1980) and 1400MHz (Stmebnng etal 
1984) The aggregate linear polarization is also Gaussian m form, occupies the central 
portion of the profile, and is never large — gomg from perhaps 35% at 430 MHz, to 25% at 
1400 MHz, and 5% at 1055 GHz. The edges of the profile are thus almost completely 
depolarized, but slightly more so on the leading than the traihng edge so that the residui 
linear ‘component’ looks a little delayed relative to the total-power profile 
A small amount of circular polarization is also observed in this pulsar’s average 
profile, usually less than 5%. The published profiles often show a little left-hand 
(positive) circular on the leading edge, changing to a little more right-hand (negative) on 
the trailing. However, other unpublished observations suggest considerable variation 
around this general configuration. 
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With not only an MP, but an IP and PC as well, the linear polarization angle (PA) 
can be measured over a large part of the star’s rotation cycle. Backer & Rankin (1980), 
however, found the linear PA behaviour of this pulsar to be unusually complex, 
apparently requiring two different angle systems. Using the information in Rankin 
& Benson (and correcting for ‘90° flips’ at about — 8 and + 12°), Lyne & Manchester 
(1988) have obtained a value of 80° for a, the angle between the magnetic and rotation 
axes of the star, by fitting the single-vector model (Radhaknshnan & Cooke 1969, 
Komesaroff 1970) to the PA data Better resolved observations are to be found in 
Rankin et al (1989) and Blaskiewicz et al (1991), and the rapid traverse of the PA near 
the peak of the profile can be seen clearly The latter study (also correcting for a ‘90° flip’ 
at about 1 70° in the 430-MHz observation) obtained values of 79° ± 1 ° and 9 1 ° ± 600° 
at 430 and 1418 MHz, respectively 

Fluctuation-spectral analysis yields a ‘weak, broad feature in the range of 0 05 to 
0 23 cl Pi with a constant strength of 0 35 . The measured center frequency is 

0 \4c/Pi (Backer 1973). In other studies (Taylor & Huguenin 1969; Lang 1969), the 
feature appeared stronger and at a frequency of 0 18 to 0.20c/Pi We note here [and 
Rankin (1986; hereafter Paper III) did not'] that there is considerable fluctuation 
power at frequencies below 0 lOc/Pi , suggesting in the terms of Paper III both core and 
conal-associated fluctuations. Finally, there is no published evidence for null pulses, 
although the limit of < 5% (Ritchmgs 1976) is not a very strong one. 

As regards classification, pulsar 0823 -1- 26 has had a chequered history In Paper III 
It was seen as a conal double (D) pulsar in an attempt to make sense of the MP and PC 
as a conal component pair The spacing between the two components was uncomfor- 
tably large, and there was some antisymmetric MP circular polarization, but the lack of 
MP conal ‘outriders’ at high frequency tilted toward the conal side The growing weight 
of evidence regarding core-component geometry in Paper IV shifted 0823 + 26 into the 
core-single (S,) class. Its MP width was exactly what would be expected for a core 
component in a pulsar with a nearly orthogonal magnetic geometry, a geometry which 
seemed highly compatible with its observed interpulse. Its rather high also 

seemed to suggest that it is not a cone-dominated pulsar And viewing the MP as 
pnmarily a core component was more compatible with its polarization characteristics 
This classification left only two nagging questions a) In a pulsar which can be observed 
to such high frequencies, where are the high frequency conal ‘outHders’*^ And b) If the 
MP IS a core component, how are we now to interpret the PC? 

These questions are the starting point for this paper In the next section we introduce 
the observations considered m our analysis, and §3 begins our discussion of the 
polarization structure of pulsar 0823 + 26’s profile. In § 4 we consider mode-separation 
techniques and their application to pulsar 0823 + 26, § 5-§ 7 explore the pulsar’s 
main-pulse, postcursor, and interpulse emission geometry, respectively, § 8 then returns 
to the polanzation-modal structure of the main-pulse emission, §9 assesses the 
significance of the displaced modal emission, and § 10 provides a sununary and the 
conclusions of the paper. 


2. Observations 


The single-pulse observations used in our analysis below come from three programs 
carried out at the Arecibo Observatory over a long period of time. The older 430-MHz 
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observations were earned out in January 1974 with a single-channel adding 
polarimeter of 500-kHz bandwidth and l.O-ms integration time, giving a nominal time 
resolution of about 0 T longitude. The polanmetry scheme is desenbed in Rankin et al 
(1975); no attempt was made to correct the Stokes parameters for the known 0 25% 
cross-coupling in the feed, which can produce spunous circular polanzation at 
a nominal level of about 10% of the linear 
The older 1400-MHz observations were carried out in October 1981, again with 
a single-channel, adding polanmeter Here the bandwidth was 10 MHz and time 
constant 0 66 ms, giving an effective resolution of about 0.5° A serious effort was made 
for the first time to correct the measured Stokes parameters for mstrumental distortion 
using the ‘orthogonal’ approximation described by Stinebnng et al (1984) 

The newer observations at 430 and 1404 MHz were made in a single observing 
session in October 1992, the 430 MHz on the morning of the 20th and the 21 cm on the 
morning of the 15th Both used a special program to gate the 40-MHz correlator, and 
the basic data recorded at the telescope were the ACFs and CCFs of the right- and 
left-hand channel voltages. The higher frequency observations used a 20-MHz B W and 
the lower a 10-MHz BW, and the retention of 32 lags m both cases reduced dispersion 
delay across the bandpass to negligible levels The resolution was then essentially the 
correlator dump time, which was 506 ^s or 0.34° The Arecibo 40-MHz correlator is 
desenbed by Hagan (1987) and the observing software by Perillat (1988, 1992) The 
measured correlation functions were scaled, 3-level sampling corrected, and Fourier 
transformed to produce raw Stokes parameters, which were in turn corrected (channel 
by channel) for dispersion, Faraday rotation, mstrumental delays, and all of the known 
feed imperfections as determined by full-sky tracks of pulsar 1929 -I- 10 and other 
sources. Dunng the course of our analysis, we discovered that the instrumental 
polarization is highly frequency dependent, particularly at 430 MHz, therefore, the 
most recent observations represent some of the best calibrated ever made at the 
Arecibo Observatory The instrumental details will be described m a forthcoming 
paper (Rankin, Rathnasree & Xilouns 1996) 


3. The polarization structure of pulsar 0823 + 26’s profile 

Figures la-c give a collection of profile and polanzation-angle density plots for both 
430 and 1404 MHz Notice the difference m the longitude scales, at 430 MHz the entire 
MP-PC portion of the rotation cycle of the pulsar is included, whereas at 1404 MHz 
only the MP region is plotted We should say that the lower panels are not really 
grey-scale plots; each point represents the polarization of an actual measured sample 
falling above a threshold of a few standard deviations in the off-pulse noise. The regions 
of greatest pnmary-mode power are thus ‘over exposed’ in the process of plotting 
enough individual pulses to explore the overall space represented by the diagram. 

These plots are interesting and remarkably complex. Let us first consider the simpler 
1404-MHz diagram (Fig. Ic), wherein it is clear that there are two distinct modes and 
that these have very different PA traverses One can also see this, both in Stinebnng 
et fl/.’s Fig 5 and in Gil et al ’s Fig 2 - we are dealing here with the same data - but 
considerably less clearly owing to the details of plotting or analysis. The darker pnmary 
mode has a slope of 14-15°/°, and the less prominent secondary mode a slope just about 
twice this What is most stnkmg is that the two modes are anything but orthogonal, and 
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near the trailing edge of the profile (near + 3° longitude) appear to be separated by 
barely 45°. Lest anyone wish to attribute this to instrumental causes, we have another, 
shorter 1404-MHz sequence from our 1992 October observations which shows an 
identical behaviour Other than this it is interesting to again note the minor peak in the 
linear power as well as the larger one, which trails the profile center Also one can see, 
both how narrow the longitude interval over which the PA can be traced is and how 
poorly It often traces the PA of the pnmary-mode power. 

We see a good deal more in the 430-MHz figures (Figs, la and b), and it is 
noteworthy how much more similar the PA density plots are than the average 
position-angle behaviour. The primary and secondary modes are again prominent in 
the center of the profile, and their measured slopes are nearly the same as at 21 cm. 



LONGITUDE [degrees] 


Figure 1(a). 
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Figure 1(b). 

Again there is rather little aggregate linearly polarized power outside the MP center, 
but owing to the greater intensity of the pulsar at 430 MHz, individual samples with 
sufficient power to define the PA can be found over nearly a quarter of the rotation 
cycle of the pulsar The principal primary-mode ‘blotch’ begins at about — 12° 
longitude and ± 90° PA and rotates positively through about 90°, also making the 
largest contribution to the overall hnearly polarized power, L Several other ‘patches’ of 
significant linear polarization can also be seen Note the area associated with the 
postcursor, another on the extreme leading edge of the profile (which can also be seen in 
L), and a third, secondary-mode ‘blotch’ which begins at about — 14° longitude and 
— 45° PA, again rotating positively through 90° At this point this last traverse appears 
to flatten off adjacent to the principle pnmary-mode ‘blotch’, and its power is conflated 
with the mam primary-mode ‘blotch’ in L. 
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Figure 1(c). 

Figure 1 (a-c). Average profiles and polanzation-angle density plots for pulsar 0823 + 26 at 
430 MHz (a, b), and 1404 MHz (c). The upper panel plots Stokes parameter /, the total power 
(solid curve), L = (Q^ + the total linear polarization (dashed curve), and V, the circular 

polanzation (dotted curve) The lower panel gives the polanzation position angle of each sample 
falbng above some threshold (see text), and the average position angle is also plotted, where 
defined, as a sohd curve The 430-MHz observations were recorded on 6th January 1974 and 
20th October 1992 and represent 3194 and 1800 pulses, respectively. The 1404-MHz observation 
was made on 9th October 1982 and represents 6200 pulses 

Again It IS interesting to consider how different the overall linear PA traverse is m the 
two 430-MHz observations, given the striking similanty of their modal-power distribu- 
tions. It does appear that the polarization of this pulsar shows significant time 
variability; the differences in fractional linear polarization, PA traverse, and the entirely 
different form of the circular are quite typical of other (not shown) observations in our 
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P0823+26 




Figure 2. (a) Average profile of pulsar 0823 + 26 at 430 MHz as m Fig 1 b At longitudes largei 
than — 15° all the quantities/, Land K have been multiplied by a factor of 5 to show the region o 
the postcursor more clearly, (b) The fractional linear polarization, <L//> (solid curve) anc 
<^>/<-f> (dashed curve). In both the curves the averages have been corrected for statistica 
polarization In regions of low intensity, the degree of polarization could be quite undefined and 
oscillating between high and low values as m the left hand portion of the figure 
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P0823+26 



Figure 3. Pulse intensity histogram (solid line) for pulsar 0823 + 26 at 430 MHz corresponding 
to the October 1992 observations in Fig lb. The dashed line shows the mtensity histogram 
obtained from the noise window, which is normalized to 142 pulses. 
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possession. It is quite conceivable that this pulsar exhibits some yet unrecognize 
mode-switching phenomenon (and of course this would be particularly difficult 1 
identify in a pulsar with a smgle profile). 

Turning now to Fig 2, we agam see the 430-MHz profile although displayed 
a slightly different manner than m Fig. lb In the top panel, at all longitudes later th£ 
~ 15°, /, L and V have been multiplied by a factor of 5, so that the posLcursor ar 
interpulse windows can be seen more clearly. The information in the bottom panel hi 
been obtained after rejecting sample by sample all data points with I < twice the a 
the noise wmdow Here, we plot the fractional linear polarization, first in the usu 
manner by averaging Stokes parameters Q and U (dashed curve) and, second, \ 
aggregating the hnearly polarized power irrespective of the PA (solid curve). Under tl 
MP peak the latter is larger than the former as expected, but generally only by son 
10% or so, arguing that most of the depolarization is occurrmg on time scales sho 
compared to the 50Q-fis resolution of the observation However, the amount 
depolarization takmg place over tune-scales larger than a penod (which amounts to tl 
diflTerence between the two curves in the bottom panel) is considerably greater in tl 
wings of the MP profile This suggests that different processes and tune scales 
depolarization are operative m the MP wings, and overall, that the MP has a composi 
character. We will return to this issue below. Finally, note that exactly the opposite 
true for the PC; the two linear polarization measures are identical within the non 
indicating that virtually no depolarization is occurnng on time scales shorter than t] 
sampling interval 

Finally, Fig. 3 gives a 430-MHz, pulse-mtensity histogram, scaled to the averag 
profile intensity to show the range of fluctuation (solid Ime) A similar histogra 
(dashed hue) for the noise window has been normalized so that the zero-intensity bi 
coincide Interestingly, the pulse-intensity histogram is bi-modal, and the distin 
zero-energy peak has a width almost exactly equal to that of the noise distnbution T 
fraction of null pulses from this analysis is 6 4 + 0 8 % as compared with the upper Im 
of 5% from Ritchmgs’ (1976) study 


4. Polarization-mode separation 

Several efforts have been made to separate the polarization modes and to reconstri. 
the two modal pulse profiles Cordes et al. (1978) first carried out this analysis for PS 
2020 + 28, and Rankin etal (1988) and Rankin (1988) then applied it to pulsa 
1737+ 13 and 1604-00, respectively. About this time Gil etal (1991) develop 
a different polanzation-mode fitting technique and applied it first to 0823 + 
observations at 1404 MHz and then to observations of a small group of other pulsars 
the same frequency (Gil et al 1992). The first techmque rehes on a model of t 
polarization angle as a function of longitude, usually the ‘single- vector’ model given 
Radhaknshnan & Cooke (1969) and extended by Komesaroff (1970), whereas t 
second uses an algonthm which partitions the aggregate polarization angles into t\ 
groups statistically, longitude by longitude. 

The two techniques are complementary in that, both have their strengths ai 
weaknesses The first technique is more sensitive and ‘stable’, but risks giving the fa^ 
appearance that the modal PAs follow the model in regions of longitude where t 
modal power is negligible. Gil’s technique, on the other hand, has the advantage 
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being model-independent, but requires significant power in both modes at each 
longitude and easily loses track of the identity of the modes when they are discontinu- 
ous Our study of 0823 -h 26 uses only the first technique, because we were unable to 
obtain sensible results using Gil’s technique on the 430-MHz observations. 

The question then immediately arises regarding how this overall polarization PA 
behaviour should be interpreted At 430 MHz, beginning with the primary mode, we 
see a positive, about 1 20° rotation of the PA which seems to terminate at around the PA 
of the PC; by this interpretation, however, there is no appreciable pnmary-mode power 
at longitudes earlier than that of the primary-mode ‘blob’ — a circumstance which is 
also largely true, we note, at 21cm What power there is at earlier longitudes is 
reasonably associated with the secondary mode, although the leading-edge patch 
seems somewhat disconnected from the mam area where the secondary mode rotates 
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Figure 4(a). 
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Figure 4(b). 

most rapidly. It furthermore seems that the primary and secondary mode 
barely distinguishable in the ‘patch’ near - 5 to - 10° longitude and 
Conveniently - despite the uncorrected ionospheric Faraday rotation - 
MHz observations have nearly the same PA origin. 

Using this interpretation of the polanzation-mode behaviour, we cai 
separate the composite profiles into their constituent polarization-mo< 
profiles. In Rankin et al (1988) we found a convenient representation of 
vector model 


- (^o) 


A-[A-j ]cos{<P-^o) 


X = Zo + tan 
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Figure 4(c). 

Figure 4(a-c). Partial, ‘mode-separated’ profiles for pulsar 0823 + 26 at 430-MHz on 6th 
January 1974, (a) pnmary-mode, (b) secondary-mode, and (c) residual profiles. The lower panels 
in these figures give the PA traverses correspondmg to the two modes, and the model PA is 
overplotted as a dashed curve 


dx 


where x is the polarization angle, (j) the longitude, R = 1^1 = sin a/sin P the maxi- 

ld0|mw 

mum PA rate, and A a constant which can be computed from the total PA traverse Xc 
between reference longitudes ± as 

sin 0c Aan^Xc-cos(^c/|K| 

1 — cos <f>i. 


A=- 


and which is formally equal to sm C/tan a, where a and C are the angles that the magnetic 
axis and the sight hne make with the rotation axis, respectively 
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Now taking the center of the prun ary-mode PA traverse at about ~ 4° 1- 
a PA of — 40°, R as about -j- 15°/°, and computing A from a total traverse 
± 30° longitude, we proceed to construct partial (‘mode-separated’) pJ 
spondmg to the primary and secondary modes of 0823 + 26 The stream « 
the individual pulses is compared with a threshold - the criteria being a) 
total linear polarization exceeds twice the off-pulse baseline noise level, an 
the computed PA falls within ± 45° of the model PA - and accumulated ir 
partial profile ‘bins’ accordingly The first, ‘pnmary-mode’ partial profih 
those samples which both exceeded the noise threshold and had PAs lyin^ 
of the Poincare sphere closest to the Stokes vector representing the me 
second those in the other hemisphere, and the third those which did not m 
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Figure 5(a). 





Polarization- Mode Separation and Pulsar 0823 + 26 


341 



-45 -30 -15 0 15 30 45 

LONGITUDE [degrees] 

Figure 5(b). 

threshold. The ‘mode-separated’ profiles at 430 MHz computed from the earlier 1974 
and the later 1992 observations are given m Figs. 4 and 5, respectively. 

Turning first to the earlier observations in Fig 4, we see that the pnmary-mode 
partial profile (Fig 4a) greatly resembles the total profile in Fig. la. The linear 
polarization is only a little greater, the circular polarization almost identical, only we 
now have an MP which is slightly narrower with its linear polarization more symmetri- 
cally positioned in its center. With a noise threshold of about l-a* only about 9.9% of 

^The precise off-pulse noise level is difficult to determine in the older observations as the 
noise-window data is no longer readily available, comparison with the 1992 observations - 
where the noise levels were accurately determined - however, shows that our estimate of c is not 
far in error 
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Figure 5(c). 


Figure 5(a-c). Partial, ‘mode-separated’ profiles for pulsar 0823 -f 26 at 430-MHz on 
October 1992 as in Fig 4 

the samples accrue to the pnmary-mode partial profile, but these carry nearly 67.9‘ 
the total power The only big change is in the PA traverse, which now more ne 
conforms to what is expected on the basis of the smgle- vector model. Part of this i 
artefact of the ‘mode-separation’ procedure; in regions of longitude where prim 
mode L is small compared to the noise level - such as the ‘baseline’ region on 
MP - the precise value of the model PA does affect the PA of the modal partial prc 
However, m regions where any significant linear power is present, the partial pn 
PAs are only very weakly dependent upon the PA model We also note that the P 
most promment in the primary mode, and its linear polarization appears ir 
complete because weakly polarized samples have not reached the noise threshold le 





Polarization-Mode Separation and Pulsar 0823 + 26 


343 


It IS the secondary-mode profile in Fig. 4b which holds the surpnses, though it 
consists of only 4.1 % of the samples and 17 7% of the power. This profile shows a more 
complex structure, both in L and in total power. Were it a complete profile, we would 
immediately say it had a tnple form The polanzation is low, but not insignificant, 
throughout its entire extent, here we see much less severely depolarized ‘wings’ as 
compared with the primary-mode profile Part of this is again an artefact of the 
separation in that, weakly polarized samples in the wings will have a much reduced 
chance of meeting the threshold in L; however, as we shall see below, this is not the 
whole story. Again, the PA traverse is ‘rationalized’, but it is also well defined over most 
of its extent by ‘patches’ of modal power. Note that here we see a flattening of the steep 
traverse after about — 15° longitude as compared to the density plot in Fig. 1 because 
the PA model cannot track both the primary- and secondary mode traverses in this 
region simultaneously 

Returning now to the form of the secondary-mode profile, let us begin to ask what 
interpretation should be made of it The first thing to point out is that the pnmary- and 
secondary-mode profiles are not coincident in time The two partial profiles have 
virtually identical half widths (fwhm) at 4.1°, however, the secondary-mode one arrives 
earlier by about 1.5°, Then there is the question of the ‘outriders’! The leading feature 
can be traced to the isolated ‘patch’ at about — 23° longitude and — 30° PA in Fig. lb 
which we pointed out earlier. The trailing feature is more difficult to see as it is below 
and to the right of - but immediately adjacent to - the principle primary-mode ‘patch’ 
at — 9 to — 14° longitude and ^ 90° PA. Note the asymmetry of the pnmary-mode 
feature, and the additional emission to the bottom right of it Note also that the 
minimum between the ‘MP’ and the trailing feature in the secondary-mode partial 
profile falls just after — 7.5° longitude — and that there is a slight ‘cusp’ in the density 
distribution in Fig. lb which apparently corresponds to it Here the PA model used in 
the mode separation is sensitive to the position of the modal boundary, so that it 
effectively discnminates between the above primary-mode ‘patch’ and the secondary- 
mode power immediately adjacent to it. This criticality stems from the close separation 
in angle of the two modes in this region 

Readers may also question whether the secondary-mode features are in some 
manner generated by the mode-separation procedure After all, most of the samples 
(86.0%) and a significant amount of the power (14.4%) in the onginal pulse sequence 
are not represented in Figs 4(a and b) and accrue to the residual profile in Fig. 4c In 
order to explore this possibility, partial profiles were computed holding the model PA 
parameters constant, but progressively reducing the noise threshold from lo to Ic, Oc, 
and even — la. Eventually, with the latter, the noise threshold becomes moot and all 
samples accrue to either one or the other modal partial profiles. The meamngless 
division of ever more noise-dominated samples into one or the other ‘modal profiles’ 
produces statistical travesties (total linear polanzation L > /), but the form, bound- 
aries, and relative amplitudes of the secondary-mode ‘MP’ and its ‘outriders’ remain 
intact. Therefore, we cannot trace their ongin to an artefact of the partial-profile 
construction process If there is an error, it must be in the interpretation of the pnmary- 
and secondary-mode related ‘patches’ in Fig. 1. 

Turning now to the more recent 1992 observations, we see a very similar story overall. 
The pulsar appears to have been about 4 times weaker on the latter occasion, but 5 times 
the bandwidth and much better receivers were also being used in 1992. It is worth 
returmng to Fig. 1 and again carefully comparing the two 430-MHz observations. The 
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MP in the 1992 observation is significantly narrower than that in the earlier one (3.7® 
versus 4 5®), owing perhaps to improved resolution. This seems to accentuate the 
asymmetry of L within the overall I envelope, and probably accounts in part for its 
larger relative amplitude as well. Very different, however, is the circular polarization; 
while the 1974 observation shows a small amount ( < 5%) which varies antisymmetn- 
cally across the MP, the 1992 profile has nearly 15% RCP with only a hint of a sense 
reversal. 

The pnmary-mode, secondary-mode, and residual partial profiles are given m 
Figs. 5(a-c), respectively Little needs be said about the former except to note its 
narrowness (now about 3.2° fwhm), its increased linear symmetry as compared with the 
total profile, its nearly identical circular polarization, and of course its more orderly PA 
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traverse With a l-a threshold, it here represents 10 1 % of the samples and 68.9% of the 
total power in the individual-pulse sequence. 

The secondary-mode profile in Fig. 5(b) again has a tnpartite form, both in 1 
and in total linear, L — although the mimmum between the middle and trailing 
features is less well resolved The fractional linear polarization seems to be a little 
larger than in the earlier observation. However, please note that the secondary-mode 
circular is very comparable, notwithstanding the very different pnmary-mode 
circular (again suggesting that a ‘mode-changing’ phenomenon is active in the pulsar 
which changes, over time, the relative proportions of the two polarization modes 
or their overall polarization states) The fwhm of the ‘MP’ is here about 3.3°, and it 
again leads the pnmary-mode one by almost exactly 1.5°. Here the secondary-mode 
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Figure 6(c). 

Figure 6 (a-c). Partial, ‘mode-separated’ profiles for pulsar 0823 -I- 26 at 1404-MHz on 9th 
October 1981 as in Fig 4 

profile represents 4 0% of the samples and 19 0% of the power of the total individual- 
pulse sequence at an L-threshold value of 2a. And, again, lowering the L threshold 
to irrelevance produces no qualitative change in the form of the secondary-mode 
profile 

‘Mode-separated’ profiles corresponding to the 1404-MHz observation in Fig. 1(c) 
are given in Figs. 6(a-c). Only a limited longitude range around the MP was available 
m these observations. The PA-model parameters used were similar to those given 
above. At this frequency 60 1% (16 5®/o) and 24.1% (8 4%) of the power (samples) 
accrued to the pnmary (Fig 6a) and secondary (Fig 6b) modes, respectively, and 
15.8% (75.1 ®^) to the residual (Fig. 6c) profile 
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Rather few ‘mode-separated’ partial profiles, computed from well resolved and calibrated 
observations, are available m the hterature; it is thus far too early to attempt any generaliza- 
tions regardmg them. In particular, the question of what constitutes a ‘component’ m total- 
average profiles has been a subject fraught with controversy (and is perhaps only now 
beginning to clarify); therefore we certainly have no wish to hastily attnbute this status 
to every discernible minor feature in mode-separated partial profiles. Nevertheless, the 
appearance of the additional features m pulsar 0823 + 26’s secondary-mode profile is 
so extraordinary that we are emboldened to attempt some speculative analysis. 

Therefore, let us proceed by entertaining the hypothesis that the two new, minor 
features on the leading and trailing edges of the mam one are indeed correctly 
interpreted as a conal-component pair, which was too weak to detect in the total profile 
at any frequency We noted above — and reiterate here - that 0823 -h 26, though almost 
certainly correctly classified as a core-single (SJ pulsar, produces no identifiable conal 
‘outnders’ even at the highest frequencies of observation. Thus, were we to interpret the 
secondary-mode features as ‘outnders’, we would encounter no conflict with our 
existing knowledge of this star 

A way to assess whether this interpretation makes sense is that of exploring its 
geometncal consequences. The basic magnetic geometry of this pulsar seems well 
established. Paper IV argued that core-component widths can be used to determine the 
latitude of the magnetic axis a, and pulsar 0823 -h 26 was discussed not only as an 
example of a pulsar whose interpolated 1-GHz core width implies that a ^ 90°, but one 
whose overall PA traverse supports this, and whose mterpulse (0823 -h 26 is almost 
certainly a two-pole interpulsar) leaves virtually no doubt. 

More questionable then is its emission geometry in relation to our sightline. 
Referring again to Fig. 1, recall that the maximum PA sweep rate R is about -j- 15°/° for 
the pnmary mode and about -h 30°/° for the secondary one. Using the definition of 
R followmg equation (1), we see that this results in P values of -1-3.8° and -h 1.9°, 
respectively, where the ‘H-’ sign has been added to mdicate that these are outside (positive) 
traverses on the basis of their flattened wings (see Narayan & Vivekanand 1982). 

The significance of the P value associated with the primary mode is not yet clear, but 
we are immediately at liberty to use the secondary-mode p value, along with a and the 
(outside, half-power) width of the conal component pair to calculate the conal emission 
radius which would be expected. This latter value, scaled from the older and newer 
secondary-mode profiles in Figs 4 and 5 are 16.7° and 16.3°, respectively. Now using 
equation (4) of Rankin (1993a; hereafter Paper Via), this implies a conal radius p of just 
less than 8 5° (if we ‘round’ down to 16° in an effort to ‘extrapolate’ to 1 GHz, p is 8 2°; 
please see Rankin (1993b); hereafter Paper VIb, Table 4). In either case, 0823 -h 26’s 
secondary-mode p value is clearly just what one would expect for an outer cone, as can 
be seen by refernng to equation (5) of Paper Via 

It would then seem that the most economical interpretation of the ‘outriders’ in the 
secondary-mode profile of pulsar 0823 + 26 is that they indeed represent an outside 
conal component pair, which is so weak that it can be detected only by removing the 
dominant primary-mode power from the profile. This line of interpretation is entirely 
consistent quantitatively with the known geometry of the pulsar and further an answer 
to the question - for 0823 -h 26, certainly, but also in general - of why some core-single 
pulsars have no discernible conal outnders 
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6 . Posteursor emission geometry 

First, several different lines of argument now give similar conclusions regarc 
basic geometry of pulsar 0823 + 26 Its interpulse (Hankins & Fowler 1986; Pa 
analysis of its PA traverse (Lyne & Manchester 1988; Blaskiewicz et al. 1991), . 
width of Its MP core component (Paper IV) all indicate that its magnetic axis is 
orthogonal to its rotational axis Moreover, two somewhat different analyse! 

& Manchester (see especially, their table 2) and that m both this paper and Pap 
find that our sight line passes within a mere 1 or 2° of the magnetic axis 

Second, the identification above of the secondary-mode ‘outnders’ as an oute 
component pair marks these components as falling at the extreme boundary o 
emission in this pulsar, in terms both of angle from the magnetic axis and of di 
from the star We now have overwhelming evidence that conal emission is 
consistent in its angular beaming properties (Paper Via & b; Gil et al 1 993, Gil & 
1993), the outer conal emission radius is 5 75° or, as we saw in the pr< 
section, about 8° 

In this geometric context it is now much more interesting to consider the signifi 
of the postcursor component in pulsar 0823 -f 26, and the first point to make is tl 
now have no possible option of viewing the PC component as some kind of 
widely separated portion of a conal double profile as was attempted in Paper II 
PC lies much too far from the magnetic axis - some 30° - for it to be of conal ori 
any usual sense; furthermore, there is only a PC, no precursor^ 

While we cannot say how it is that the PC is emitted, it would seem that 0823 
geometry offers a particularly good opportunity to identify unusual types of emi 
As noted above, not only is this pulsar’s magnetic moment nearly parallel to its eqi 
but so is our sight line' Only in this geometry is there a bundle of field lines (also p£ 
to the star’s rotational equator) which remains in view as the star rotates an 
tangent points in the observer’s direction for an unusually large fraction of the 
rotation period One such bundle trails the magnetic axis, and it is worth noting 
any ‘field-line sweepback’ at high altitude will increase the curvature of this bi 
Another bundle leads the magnetic axis, and (with enough sweepback) it coulc 
have tangents in the observer’s direction, but only at great altitude and for a sh 
interval 

If the PC component is emitted along this trailing bundle of peripheral, equal 
field lines at a height at which the field remains essentially dipolar, then it is possil 
estimate the emission height of the PC, again using equation (2) in Paper Via F. 
angular separation between the PC and the magnetic axis (taken to be the MP pet 
about 32° (Hankins & Fowler 1986), we find an emission height of about 4 x 10^ 
that is, about 400 stellar radii or some 4% of the light-cylinder radius - as comi 
to the nominal 210-km emission height of the outer cone. Any poloidal distortK 
the field would increase the curvature of the trailing bundle of equatorial field 
and decrease the emission height Interestingly, it would also decrease that o 
leading bundle, perhaps explaining why we more generally see postcursors n 
than precursors. Both height and sweepback will tend to increase the curvatu 
the trailing bundle, which may favour emission by particles with quite diffe 
energies than those responsible for the MP Theories of pair-production disch 
above pulsar polar caps, seem to produce a broad distribution of particle enej 
(Daugherty and Harding 1986) Finally, we note that the MP-PC separation mere 
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weakly with frequency, suggesting that the higher frequency PC emission is emitted at 
greater heights than the low frequency emission, if so, this would be a ‘radius-to- 
frequency mapping’ of the opposite sense to that usually encountered 

If this is the correct interpretation for the geometry of the PC component, then we 
note that both the PC and the latter conal outrider are emitted along the same bundle 
of trailing equatorial field lines It is then possible that the two components will exhibit 
some correlation in the intensities of their constituent individual pulses The longitude- 
to-longitude correlation of the total power, correspondmg to the 430-MHz observa- 
tions in Fig lb, IS given in Fig 7. The longitude, of course, labels both axes, and average 
profiles composing the noise, the MP-PC, and IP windows have been plotted for 
convenience. A colour scale indicating the degree of correlation is given at the upper left 
of the figure The upper tnangular region shows that the longitude-longitude correla- 
tions are calculated by correlating the emission in the same pulse at different longitudes 
(zero lag) In the lower triangular region, the emission at any given longitude corre- 
spondmg to the row number is correlated with emission in the previous pulse at 
longitude corresponding to the column number (lag= + 1) The diagonal line is 
included m the lower part with lag = + 1 

This figure clearly shows the correlation between the MP and PC emission 
Interestingly, the correlation is significant only between the central portion of the PC 
and the wings of the MP, between the center of the MP and the PC there is hardly any 
significant correlation. If anything, the emission m the center of the MP is anti- 
correlated with the PC emission of the previous pulse, as can be seen from the lower 
part of the figure However, the correlation between the wings of the MP and PC 
persists even for a one-penod lag Similarly, the MP leading and trailing wings exhibit 
highly correlated emission, but this emission is not significantly correlated with that of 
the MP center Another curious feature is the anti-correlation between the MP peak 
and the ‘bridge’ region between the MP and PC. This correlation was also noted by 
Popov & Sieber (1990) in their 1700-MHz study, although their observations were 
restncted narrowly to the MP region There is also a significant correlation between the 
leading part of the previous MP and the current IP peak 

AU this notwithstanding, we note that there are five other pulsars with known 
PC components: The Crab pulsar is usually discussed as having a ‘precursor’ compo- 
nent, but in that this component appears to be a core component (Paper IV), one might 
better regard the star’s MP as a postcursor to it Then there are 0940 + 16, 1530 + 27 
(Rankin etal 1989), probably 1530 — 53 (McCulloch etal 1978, Manchester etal 
1980; McCulloch etal 1982), 1929 + 10 (Perry & Lyne 1985, Rankin & Rathnasree 
1996), and more recently 2217 + 47 (Suleymanova & Shitov 1994). We know of no 
pulsar with a true precursor - that is, one with a weak patch of emission which 
leads a core component Apart from the Crab pulsar and 1929 + 10- which are 
almost certainly orthogonal rotators - we know very little about the geometry of 
these other stars; most have not even been adequately classified, in part because of 
their odd ‘postcursor’ emission Somewhat more is known about 2217 + 47, it was 
classified as an pulsar in Paper VIb with an a of only 42° This value could easily 
increase, however, as the profiles available in the literature are mostly poor (the only 
-^400-MHz profile available is from LSG m 1971^) and thus probably not well 
resolved 

We note that pulsars 0823 + 26 and 2217 + 47 have very similar periods and period 
derivatives, and the completely equatorial model used above to discuss the emission 
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geometry of the former also provides an attractive scenario for understanding 
Suleymanova& Shitov’s observations of the latter. Perhaps the most surprising feature 
of the 2217 -h 47 postcursor is that it slowly varied in amplitude and position over time, 
possibly correlated with changes in the torque on the star Indeed, if the PC emission is 
occurnng high on this bundle of trailing, equatorial field lines, torque changes might 
change the curvature and thus the height of the radiation, which would in turn change 
Its phase within the profile relative to the core component. But we do not yet know 
whether 2217 + 47 -or really any of the other ‘postcursor’ pulsars apart from the 
Crab -have orthogonal geometries We are guessing that they probably do, but 
additional observations and analysis are required to answer this question If so, 
some - 2217 + 47 to be sure - may prove highly useful as a means of exploring other 
regions of the pulsar magnetosphere. Pulsars with PC components may also be good 
candidates in which to search for new interpulses 


7. Interpulse emission geometry 

Little can be said with confidence about pulsar 0823 + 26’s interpulse geometry 
No quality observation of this feature exists in the literature, and, given its weakness, 
none may be possible in the foreseeable future. Nonetheless, the 430-MHz profile 
(Rankin & Benson 1981) appears more single, and the 1400-MHz profile (Rankin et al 
1989) has a flat top which might indicate incipient conal outriders (as per the usual 
tnple or possibly conal-single evolution) Both IP profiles have a full width at half 
maximum of about 16 or 17°, and the PA sweep rate is about — 67°. These values do 
not make any immediate sense in terms of the geometrical arguments of Paper Via & b 
If the above width corresponds to an outer cone, either the IP sweep rate needs to be 
steeper or the width smaller One possibility is that here we are cutting an outer cone on 
Its extreme edge, so that the sightlme comes nowhere near the high-mtensity ‘shoulder’ 
of the conal beam. 


8. Polarization-modal structure of PSR 0823 + 26’s emission 

We have seen above that the pnmary mode at 430 MHz represents about 68-69% and 
the secondary-mode about 18-19% of the total power in the sequence of pulses when 
the threshold is set at a level such that L must exceed about two a in the ofl-pulse noise; 
that is, the primary mode is about 3.5 times stronger than the secondary one. Of course, 
these numbers will change somewhat if a) the mode-separation algorithm is not 
accurately segregating the modal power in the PA-longitude diagram, b) if the 
L threshold level is changed, or c) if the emitted modal power varies m time We have 
done our best to assure that the modal separation is correct, although given the 
proximity of the two modes near the primary-mode ‘blotch’ in Fig. 1, some ambiguity is 
unavoidable. As for the second condition, a 2-a threshold gives a sample uncertainty in 
the PA of + 30°, and any lower threshold will increase the number of unpolarized or 
noise-dominated samples which are spuriously assigned to one mode or the other We 
re-emphasize that only for those samples with adequate residual linear polanzation can 
we make a choice about their mode For the others, all such information has been lost 
through depolarization on time scales short compared to the sampling interval, and we 
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can only speculate about the significance of this depolarized portion of the puls 
emission Observations suggest either that these are relatively stable or that we h 
found the pulsar in the same ‘modal’ state 

We can then ask whether this is the ‘true’ power ratio of the two modes. So 
12-14% of the power fell below the threshold in L, and the average of these ‘residi 
samples shows virtually no linear polarization It might then be a) that most 
this power is associated with the secondary mode, which is only some 10% polanz 
m which case the modal ratio is about 70/30%, b) that the power is about equa 
divided between the two modes (and, indeed, the residual PA falls near the prim£ 
mode before the MP peak and near the secondary mode thereafter) in which a 
the modal ratio is about 75/25%, or c) that most of the residual power is associat 
with the primary mode in view of its greater intensity, in which case the moc 
ratios are about 80/20% At 21 cm about 60% of the power was found in the pnma 
mode to about 24% and 16% for the secondary mode and residual categories, s 
Fig 6 If the polarization-mode switching phenomenon is broad band - as, indee 
it appears to be -these latter numbers might be squared with any of the thr 
alternatives 

It is also useful to consider the fractional linear polarization, and companng the tot 
average profiles m Fig 1 with the modal partial profiles in Figs. 4-6 we see that tl 
mode separation did not dramatically increase the fractional linear polarization, eith^ 
at 430 or 1400 MHz. The earlier 430-MHz observation has a peak fractional linei 
polarization of about 35%, whereas the primary mode’s is only about 42%. Similarl 
that of the later 430-MHz run goes from about 45% to 52% and the 1400-MHz oi 
from about 28% to 36%. The later data are better resolved than the earlier at both 43 
and 1404 MHz, and this may be responsible, in part, for their somewhat larger amoun 
of aggregate linear polarization. 

In all cases the secondary-mode fractional linear is lower. We then apparently hav 
a situation in pulsar 0823 + 26 wherein very significant amounts of depolarization ar 
occurnng on time scales shorter than the inherent time resolution of the observations 
here in all cases some part of 1° longitude. This also largely explains the fact that th 
secondary mode is more depolarized than the primary one, because sccondary-mod 
samples will typically contain more primary-mode power than primary ones wil 
secondary-mode power, by simple reason of their relative intensities Only thosi 
unusual samples in which the secondary mode dominates the primary will be recog 
nized as ‘secondary-mode’ samples — and accumulate in the ‘secondary-mode’ profile 
Generally the secondary mode will just exceed the primary, and these are th( 
conditions which are most depolarizing. These factors are especially complicated ir 
0823 + 26, however, owing to the spectacular non-orthogonality of the two modes 
throughout the duration of the MP 


9. Displaced modal emission in the main pulse 

A very interesting feature of the modal partial profiles in Figs. 4-6 is the displacement 
in phase between the primary and secondary modes Reference to these figures shows 
that the secondary mode ‘MP’ always leads the primary-mode feature, at 430 MHz by 
1 5 ± 0.1° and at 1404 MHz by some 1 3 ± 0.1°. Displacement in rotational phase has 
been seen before in the components of mode-separated profiles - the first component of 
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1737+ 13, for instance, shows a dramatic offset (Rankin etal. 1988) -however, this 
effect has not, before to our knowledge, been identified in a core component'*. 

The question is what interpretation to make of this phase delay, and at this point two 
possibilities seem open The first is that the offset is a geometrical phenomenon We saw 
above that the observed offset was somewhat greater at 430 as compared to 
1400 MHz -but just so within the errors -and just over 1° of longitude is quite 
comparable to the 3 36“ diameter of the pulsar’s polar cap. Another possibility is that 
we are here seeing an effect of birefrmgent propagation in the pulsar’s inner magneto- 
sphere as suggested, for instance, by Barnard & Arons { 1986) If this is the case we might 
expect the effect to become more pronounced at lower frequencies — and furthermore 
this ‘double’ nature of the core component might be connected with such well known, 
but entirely unexplained phenomena such as ‘absorption’ and the observed broadening 
of some core components at lower frequencies (see Paper II). 


10. Summary and conclusions 

The polarization structure of the radiation which comprises pulsar 0823 + 26’s average 
profile at both 430 and 1400 MHz has been studied Two polarization modes are 
identified in the emission over a substantial range of longitudes around the MP and 
PC In some intervals they are roughly orthogonal; however, just under the MP the two 
modes are more conspicuously non-orthogonal - and thus difficult to distinguish - 
than in any other pulsar we know of We have attempted to construct partial profiles in 
which the power in the two polarization modes is largely separated and have identified 
two (mostly) continuous, non-orthogonal PA trajectories which seem to correspond to 
these modes. The primary mode in 0823 + 26 is three to perhaps five times stronger 
than the secondary mode overall and thus dominates the emission and the PA at most 
longitudes. It is then not surprising that the primary-mode MP closely resembles the 
total-profile MP 

The secondary-mode profile consists of three MP features, which appear to be 
a central core component and a surrounding pair of conal outriders. The (outside 
half-power) width of the outriders together with the secondary-mode PA traverse 
appear compatible with the interpretation that these conal features represent an outer 
emission cone. 

This technique could lead to the discovery of covert emission features in the profiles 
of other pulsars with single profiles but complicated PA histograms. In this context, it is 
interesting to speculate about the emission geometry of the PC feature, which follows 
the MP by some 30°. Having now identified a pair of much more closely spaced 
components, which appear to represent an outer cone of emission, it is clear that the PC 
cannot be a conal component in any usual sense 

Overall, 0823 + 26 seems to have a nearly equatorial geometry — that is, both the 
magnetic axis and the sight line appear to he close to the rotational equator of this 
pulsar We then note that in this situation there is a bundle of trailmg field lines which 
has a tangent in our direction for a substantial portion of the rotation cycle of the 

Pulsar 1 737 + 1 3’s core component is also misaligned in the mode-separated profiles, but given 
the greater complexity of both partial profiles, it is not at all clear what interpretation to make of 
this circumstance 
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pulsar Perhaps it is emission along this bundle - at a height of some 4 x 1( 
about 4% of the light cylinder - that the PC emission occurs. 

If the PC component and the trailing conal outrider are both emitted al 
trailing bundle of equatonal held lines, then it is possible that their emi 
correlated Indeed, we find that there is significant correlation between the PC 
longitude ranges of the trailing conal outriders, and furthermore that this cor 
becomes negligible between the PC and the core-emission region. 

Finally, it seems that the pulsar does exhibit some 64 + 08% null pulses, i 
somewhat more than the 5% limit of Rltchmgs (1976). 
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Abstract. For accretion on to neutron stars possessing weak surface 
magnetic fields and substantial rotation rates (corresponding to the secular 
instability limit), we calculate the disk and surface layer luminosities general 
relativistically using the Hartle & Thorne formalism, and illustrate these 
quantities for a set of representative neutron star equations of state We also 
discuss the related problem of the angular momentum evolution of such 
neutron stars and give a quantitative estimate for this accretion driven 
change in angular momentum Rotation always increases the disk luminos- 
ity and reduces the rate of angular momentum evolution These effects have 
relevance for observations of low-mass X-ray binaries 

Key words. Stars: accretion — stars: accretion disks — stars* neutron — 
stars rotation. 


1. Introduction 

Disk accretion on to a neutron star possessing a weak surface magnetic field (jB < 10® G) 
provides interesting X-ray emission scenarios, and is relevant for understanding X-ray 
bursters and low-mass X-ray binanes (e.g. van Paradijs 1991) Such weak-field neutron 
stars can rotate very rapidly and are also seen as millisecond pulsars (Radhaknshnan 
& Srmivasan 1982, Alpar et al 1982; Bhattacharya & van den Heuvel 1991) and may be 
relevant for quasi-penodic oscillators (QPOs) (Priedhorsky 1986; Paczynski 1987). The 
equation of state of neutron star matter as well as general relativity play essential roles in 
such a scenario. This is in contrast to the strongly magnetic 10^^ G) accreting 
neutron stars, where plasma processes dominate (e.g., Ghosh & Lamb 1991). For the 
weak-field case, the radius of the innermost stable circular orbit (r^^^b) plays a central 
role, deciding quantities of observational interest such as the disk luminosity. The 
relevance of this parameter (r^^^) was emphasized by Kluzniak & Wagoner (1985, 
hereafter KW), who pointed out that for weak-field accreting neutron stars it is 
incorrect to always make the usual assumption that the accretion disk extends very 
close to the surface of the star, and is separated from it by a thin boundary layer. Using 
Schwarzschild geometry, Syunyaev & Shakura (1986, hereafter SS) concluded that the 
boundary layer brushing the neutron star surface will be substantially more X-ray 
luminous than the extended accretion disk If the star’s radius {R) is less than the 
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boundary layer is likely to be characterized by poorly collimated tangential motion of t 
mfallmg matter and a comparatively soft emission spectrum Whether or not R excee 
(and consequently the detailed features of the accretion scenario) depends on t 
geometry of the spacetime and also the equation of state of neutron star matter 
An important aspect of disk accretion on to weak-field neutron star is the possibih 
that the neutron star will get spun up to very short rotation penods (<millisecon 
over a time of the order of hundreds of millions of years. For such rapid rates 
rotation, the relativistic effect of dragging of inertial frames in the vicinity of the neutn 
star will be important This effect will alter the trajectories of infalling particles 
compared to the non-rotational case Therefore, for a quantitative description of tl 
accretion features, one must take into account the relativistic effects of rotation upc 
the spacetime geometry. Although the possible importance of such effects was stress^ 
by KW and SS, the details were not worked out by them 
In this paper we address this question and calculate the disk and surface lay- 
luminosities incorporating the rotational effects in a general relativistic framework. V\ 
take the spun up neutron star to be rotating at a particular value, namely, the seculi 
instabihty limit so as to illustrate the maximal reasonable effects of rotation. Th 
corresponds to the late stages of accretion. We use the Hartle & Thorne (1968, hereaft< 
HT) formahsm, for our purpose, this formalism describes a rotationally perturbe 
Schwarzschild space-time. The HT formahsm is valid for strong gravitational fields bi 
only m the limit of uniform rotation with a rate that is ‘slow’ compared to the critic; 
speed for centrifugal break-up Neutron star models rotating at the secular instabilil 
limit (assuming the star to be homogeneous), relevant in the context of accretio 
induced spun up neutron stars, are within this limit (Datta & Ray 1983), so th 
approximation will usually be adequate Recently Cook et al (1994a & b) calculate 
last stable orbits of rotating neutron stars incorporating higher order rotational terrr 
that go beyond the HT approximation. However, these authors did not use their resull 
to estimate the disk luminosity and angular momentum evolution involved in th 
accretion scenano Although the Hartle & Thorne prescnption assumes rotationall 
perturbed geometry, the use of this prescnption provides a first estimate, which i 
amenable to a straightforward numerical treatment, of the luminosity values ani 
angular momentum evolution for an accreting weak-field neutron star. This can b 
quite useful m observational applications such as low-mass X-ray binanes. Ou 
calculations are done for a range of stable neutron star configurations computed usin. 
a representative sample of the proposed equations of state of neutron star matter W 
also consider the accretion dnven evolution of the angular momentum of the neutroi 
star in a more accurate fashion than was done by KW 


2. Accretion on to a rotating neutron star 

To describe the spacetime around a rotating neutron star we use the metric suggestec 
by Hartle & Thome (1968) This metric describes a rotationally perturbed Schwar 
zschild geometry to order where SI is the angular velocity of the star as seen b] 
a distant observer. The general form of the metric is (signature: H ) 

= 0,pdx‘dx^ (a, ^ = 0, 1,2,3) 

= e^*df^ - e^'i'(d<t> - codt)^ - e^'‘de^ - e^^dr^ + 


(1 
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Here o is the angular velocity of the cumulative dragging of inertial frames and 
the critical angular velocity for equatorial mass shedding, where 
M and R are the mass and radius of the non-rotating neutron star. For simplicity, we 
use the geometric units’ c = 1 = G The metric components correspond to an interior 
with the identification: 


e^^ = e^''{l+2iho + h^P2)}, 

(2) 

= sin^ 0 { 1 + 2(i;j - h2)P2 }, 

(3) 


(4) 

21 1 + 2(»io + - 2m) 

l-2m/r 

(5) 

(where 2v is the gravitational potential function for the non-rotating star and m is the 
gravitational mass contained within a volume of radius r) and to an exterior with the 
identification: 

M' P 

e^^ = e-^^=\-l— + 2-r, 
r r 

(6) 

e^* = r^sm^e. 

(7) 

e^>‘ = r\ 

(8) 


Here M' and J are respectively the mass and angular momentum of the rotating 
configuration of the star. The quantity P 2 is the Legendre polynomial of order 2, and 
^0 > ^2 > ^0 ’ ^2 5 ^2 functions of r that are proportional to (see HT). The metnc has 

the desuable property that the mtemal and external forms match at the surface of the 
star. For our purpose here, we shall retain only the spherical deformation terms 
(characterized by subscript 0) and neglect the quadrupole deformation terms (character- 
ized by subscnpt 2). The latter are necessary for computing stellar quadrupole deforma- 
tion, but average out in calculating the rotation induced changes to M and R The 
applicability of the metnc (1) is valid for Q small in companson to and to go beyond 
that approximation requires a treatment similar to that of Cook et al. (1994a,b). 

A relativistic effect of rotation, important for the astrophysical scenano that we 
consider here, is the dragging of inertial frames, which implies 

cD(r)?i£2, (9) 

where d>{r) is the angular velocity of the stellar fluid relative to the local inertial frame, 


and is given by (HT) 



(10) 

where 


j{r) = e“''(l — 2m/r)^^^, 

(11) 


with the boundary conditions 

= 0 ; d)(r = 00) = n. 


dd) 

dr 


r = 0 


( 12 ) 
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For r > R(i e , outside the star), 
d)(r) = £2-2y/^■^ 

where J is the angular momentum of the star. 

6 


(1 


(1 


The equations of motion are derived from the Lagrangian corresponding to the metr 

( 1 ): 


= \ - oiif }, (1 

where a dot represents a denvative with respect to proper time We shall, for purpose \ 
illustration, take the polar angle (0) to be fixed and equal to 7c/2, which corresponds 1 
the equatorial plane 

The metnc (1) and the equations of motion provide three equations in the thn 
vanables (^, t, and r. From these, we can get (see KW) 


= £2 - IcoEl -h\l+ P/r^), 

(!' 

where 


E = h^i + (or^<j), 

(1' 

and 


s 

1 

II 

(li 


stand respectively for the energy and angular momentum per unit rest mass (denoted b 
m ^) and 



h = [\-2M'lryi^. 

(li 

The conditions for the turning point of the motion, the extremum of the energy, an 
the minimum of the energy are respectively given by (see Misner, Thome & Wheek 

1974) 

E^ = V^, 

d£ „ dF 
dr dr 

(2( 

(21 

and 

d^£ « 
dr" 

(2: 

where Kis the effective potential, given by 



F" = 2cu£/ + /i"(l+P/r") 

(2- 


For marginally stable orbits we can use the condition 
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Wntten explicitly, equations (20), (21), and (24) respectively become 



2 3jaE 
^ ^ 

(2x — 3) (2x — 3) 

x^ — 3a^x H- 6(a^ —jaE) = 0, 


(25) 

(26) 
(27) 


where x = r/2M\ the dimensionless radial co-ordinate;/ = J/M'^ and a = T/2M' are the 
dimensionless angular momenta for the star and the mfallmg matter respectively. 
Solving equations (25)-(27) simultaneously gives the values for the energy and the 

specific angular momentum (a^^^) the accreted particle in the innermost stable orbit 
having a radius x = x^^b 

Using Schwarzschild geometry, it was shown by SS that an accretion disk whose 
luminosity is small compared to the Eddington limit can exist only for R greater than 
the radius of the last stable circular orbit, 3r^, where = 2M = the Schwarzschild 
radius. For R < 3r^, the accreting matter falling on to the neutron star will follow the 
trajectory of a free particle in this geometry, with an energy equal to These 

authors also gave estimates of the energy released in the disk and in the boundary layer 
for R/Vg ranging from 1 5 to 10 without reference to any specific equation of state model. 
In the case where the neutron star has a high spin rate, the innermost stable orbit was 
calculated by KW using the HT metric to describe rotating spacetime These authors, 
however, did not estimate the luminosity, and they neglected the rotational corrections 
to M and R in computing the evolution of angular momentum 

For a low-field, accreting neutron star possessing substantial rotation, the luminos- 
ity from the disk accretion can be calculated using equations (25)-(27). We use here the 
notation x* = R!/2M\ where M' and R! correspond to the neutron star mass and radius 
that include corrections due to the rotation The following distinct cases are possible: 


2 1 Case (a)- Radius of the star is greater than 

If an accretion disk were to form around a relatively large neutron star (i.e., x* > 
the ingress of a particle of rest mass from infinity to the inner disk boundary will 
release an amount of energy given by 

= (28) 

where E^{x*) stands for the specific energy of the particle in the stable orbit just above 
the surface, obtained by solving equations (25) and (26) numerically for x = x* The 
energy loss in the boundary layer wiU be 

E's = m^{E,(x*)-E^{x*)}, (29) 

where Eq{x*) is the energy of the particle at rest on the surface of the neutron star, which 
can be calculated from equation (25) for x = x* and the specific angular momentum at 
the star’s surface. 
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2 2 Case (b) Radius of the star is smaller than 

In this case, x* < and the accretion disk will extend inward to a radius correspond 
mg to X = x^^y^. Now the energy released in the disk as the particle comes in from infinit 
to the innermost stable orbit will be 

= (3C 

and the energy released in the boundary layer will be 

E's = Tn,{E„^-E^(x*)}. (31 

A remark about the boundary layer luminosity formula that we use, is in order here 
Although this formula gives a plausible estimate for the boundary layer luminositj 
strictly speaking, it is an overestimate as it does not take mto account the subtraction c 
the energy that goes into spinning up of the neutron star The need for such a correctio 
was pointed out by Kluzniak (1987) and quantitative estimates for this were suggest© 
by Ghosh, Lamb & Pethick (1977), Papaloizou & Stanley (1986), and Kley (1991 
A fairly simple and general way to estimate the same was given recently by Pophar 
& Narayan (1995) who considered the accretion disk boundary layer problem ii 
cataclysmic variables. Sigmficantly, these authors have stressed that most of th 
decrease in the boundary layer luminosity will occur in the early stages of spm-up, wha 
the star is rotating slowly, rather than the late stages when it is approaching limitin, 
break up rotation rates The scenano that we consider in this paper corresponds to th 
latter stage. Therefore, for our purpose, the formula for the boundary layer luimnosit 
that we have used is expected to be adequate This point is elaborated in §4 


3. Equation of state and rotation-induced changes in structure 

The structure of neutron stars depends sensitively on the equation of state at hig] 
densities, especially for density regions > lO^^'^gcm"^ There is no general consensu 
on the exact behaviour of the equation of state at these high densities. For our purpoa 
here we choose the following six equations of state (1) Pandharipande (N) (neutroi 
matter) model, based on the lowest order constrained vanational method using Reu 
potentials (Pandharipande 1971a), (2) Pandharipande (Y) (1971b) hyperomc matter; (3 
Bethe- Johnson model V (N) for neutron matter (Bethe & Johnson 1974) which usei 
improved phenomenological potentials; (4) Walecka (1974) model for neutron matter 
using scalar-vector interactions in a field theoretical framework; (5) Wiring£ 
UV14 + UVII model for neutron-rich matter in beta equilibnum-a variationa 
calculation mcorporating the three-body interactions (Winnga, Fiks & Fabrocm 
1988), and (6) neutron-nch matter in beta equilibrium, based on the chiral sigma mode 
(Sahu, Basu & Datta 1993) Of these, models (1) and (2) are ‘soft’ equations of state anc 
(4) and (6) are rather ‘stiff’ ones, while the models (3) and (5) are roughly intermediate ir 
‘stiffness’ 

The composite equation of state to determine the neutron star structure was 
constructed by joining the selected high density equation of state to that of Negele & 
Vautherin (1973) for the density range (10‘‘^-5 x 10^°)gcm■■^ Baym, Pethick 
& Sutherland (1971) for densities down to ~ 10^ gcm“ ^ and Feynman, Metropolis & 
Teller (1949) for densities less than 10^ gem" 
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For a fixed central density (and a chosen equation of state model), the fractional 
changes m the gravitational mass (AM/M) and radius {AR/R) of the neutron star due to 
the rotation mduced spherical deformation are proportional to (Q is the angular 
velocity of the star as seen by a distant observer), and can be (numencally) obtained from 
a knowledge of the radial distnbutions of the mass and pressure perturbation terms, (r) 
and Po(r) (HT; Datta & Ray 1983; Datta 1988). The non-rotating mass (M) and radius (R) 
are obtained by numencally integrating the relativistic equations for hydrostatic 
equilibrium (see e.g, Arnett & Bowers 1977). The changes AM and AR are given by 


AM = mo(R) + 


AR=- 


Poip + p) 


dp/dr 


|r = J? 


(32) 

(33) 


where p{r) and p{r) are the pressure and the total mass-energy density at the radial 
distance r from the center of the star. 

For purpose of estimating the angular momentum evolution of the accreting neutron 
star (discussed in the next section), we need to know the baryomc mass of the rotating 
neutron star. The rotation mduced change in the baryomc mass (denoted by AM^) is 
conveniently wntten as 


AMg = AEg + AM, 


(34) 


where AEg is the rotation induced change in the binding energy of the star (HT)- 


where 


A£b = - P/R^ + I* 4TtPB{r)dr, 

1 r 1 2 JdaiV Id/ 3.2I 
4nr^\_l2’ \dr J 3 dr J’ 


(35) 


(36) 


and e = p — rngti is the density of internal energy, with n{r) and rrig denoting respective- 
ly, the baryomc density and rest mass 

Density profiles of neutron stars are remarkably flat out to r = (0-8-0-85)i?(Datta 
et al. 1995). Therefore, the concept of rotational secular instability in the context of 
Maclaurin spheroids (Chandrasekhar 1969) is a relevant approximation when con- 
sidenng the rotational stability of neutron stars. For a uniformly rotating homogene- 
ous spheroid, this instability corresponds to an angular velocity Q = £1,, given by 


ji- 

2nGp 


018, 


(37) 


where p is the average density of the star. The quantity fi, sets a rough limit up to which 
the neutron star can be spun up, before the onset of rotational instabilities. 



Table 1. Changes in Neutron Star Properties due to Rotation . 
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The prescription mentioned earlier to calculate the mass and radius of a rotatin 
relativistic star is valid only for rotation rates that are ‘slow’ m comparison to 
Hartle & Thorne (1968) constructed ‘slowly’ rotating neutron star models all the wa 
up to Q = for the Harrison-Wheeler and Tsuruta-Cameron equations of state. Her 
we shall only consider rotating neutron stars with f 2 = because of the secula 
instability consideration, and will employ the newer equation of state models men 
tioned earlier, these models will then illustrate the maximal reasonable effects c 
rotation. Since Q, = (0*27)^^^ the models constructed by us are adequately treated u 
the hmit of ‘slow’ rotation. To treat the maximal possible effects of rotation, one has t< 
go beyond this limit, and use a techmque similar to that employed by Cook et ai (1994a, b 
for very rapidly (but still uniformly) rotating neutron stars, they determined that spin-up tc 
the millisecond pulsar regime is possible for a wide range of equations of state 


4. Results and conclusions 

In Table 1 we give the neutron star properties calculated for the vanous equations o 
state as a function of the central density (pj. Columns (3) and (4) list the non-rotating 
mass and radius Column (5) gives the corresponding value of £2^ and columns ( 6 ) anc 
(7) give the differences between the rotating and non-rotating neutron star mass and 
radius. Column ( 8 ) gives the non-rotating values of while column (9) shows how 
r^rb decreases due to rotation. Typical mcreases in mass are (5-1 1)%, in radius ^ 4% 
while r^j^b decreases between 15% and 25% 

In the second and third columns of Table 2 we give values of the non-rotating mass 
M and the rotationally enhanced mass M’ for a maximal reasonable rotation rate 
(Q = Q^) Columns (4) and (7) show the disk and boundary layer luminosities for 
non-rotating configurations (£p and respectively) The values of and £ 5 , 
including rotational effects treated consistently within the HT framework, are given in 
columns (5) and ( 8 ). All values of luminosity listed in the table are in units of the 
baryonic rest mass. The boundary layer luminosity values listed in Table 2 do not 
include corrections for the energy that goes into spinning up the neutron star. We have 
made estimate of this correction following the prescription of Popham & Narayan 
(1995) for the case of 1-4 Mq neutron star corresponding to the EOS model by Winnga 
et al (1988) This is shown in Fig. 1 which is a plot of the rotational energy correction 
^Tot of the particle rest mass vs. the stellar rotation rate £2. It is clear from this 

figure that the rotational energy correction to the boundary layer luminosity tends to 
become unimportant as the angular velocity approaches very high values. For Q = £2^, 
-^rot 0*048; this is an order of magmtude less than the boundary layer luminosity given 
in Table 2. Therefore the neglect of rotational energy correction to calculate £' for the 
neutron star rotating at the limiting £2 as considered by us does not result in a gross 
overestimate 

We find that as M' increases, the luminosity in the boundary layer increases 
monotonically. The disk luminosity increases in the rotational case as compared to the 
static case; the increase, illustrated m column ( 6 ), is most significant (> 30%) for high 
M' values. The fractional changes AE^/E^ values (column (9)), however, exhibit 
a non-monotonic behaviour with increasing M\ though they typically decrease by 
a few per cent F urther, as is well known, we see that the radius of the innermost stable 
orbit IS less than the static case (where = 6M) for all equations of state considered. 
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n 

Figure 1. Rotational energy correction E^Jmc^ to the boundary layer luminosity as a function 
of the stellar rotation rate Q. 


We now make a few remarks about the angular momentum evolution of the 
accreting neutron star In performing the angular momentum evolution calculations, 
we have made the following additional assumptions: 

■ The central density of the neutron star does not change with mass accretion. This is 
a reasonable assumption to make because for a static, stable neutron star, while the 
central density increases with increasing mass, the rotation tends to ‘decrease’ the 
central density Thus if the total accreted mass is not large, the effect of rotation will 
nearly compensate for the effect of accretion on the central density; 

■ The mass accretion rate is equal to dM^/dt as in KW, but with the difference that 
rotational effects have been taken into account, that is, = dM^/dt, 

■ The maximum angular momentum for a particular neutron star configuration is 
determined by the secular instability limit: 

Q = n^ = (0-27)'/2a,. 


KW derived the evolution equation 


dJ 

dM' 


/(r(M0) 


dM' 




and solved it under the simplifying assumptions 


M'^M, 


dM'g dMB 

dM' “ dM 


(38) 


For dMJdM, KW used a fit to the values of M and R taken from Arnett & Bowers 
(1977) Here, we solve equation (38) numencally, incorporating the effect of rotation as 
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Table 3. Gravitational and Baryomc Masses as Functions of Angular Momentum. 


EOS Model 
(1) 

pc 

(g cm"®) 
(2) 

M/Mq 

(3) 

M'/Mq 

(4) 

M'bIMq 

(5) 

j 

(g cm® s~®) 
(6) 

PandJiaripande (N) 

1 SSEIS 

1 400 

1 400 

1 567 

0 OOOEOO 




1 413 

1.578 

2.526E48 




1 426 

1 589 

3 648E48 




1 438 

1 599 

4.420E48 




1 467 

1.623 

5 823E48 

Pandharipande (Y) 

4 79E15 

1 400 

1 400 

1 615 

0 OOOEOO 




1.416 

1 627 

3 269i;48 




1.432 

1.639 

4 660E48 




1.448 

1 652 

5 772E48 




1 466 

1.665 

6 746E48 

Bethe-Johnson V (N) 

i.arEis 

1 400 

1 400 

1.556 

0 OOOEOO 




1 420 

1.572 

3.139E48 




1 440 

1 589 

4 424E48 




1 460 

1 606 

5 422E48 




1.500 

1 639 

6 992E48 

WaJecka 

8.02E14 

1400 

1 400 

1 671 

0 000£00 




1 430 

1.699 

3.760i;48 




1.486 

1.750 

6 318E48 




1.500 

1.764 

6.844E48 




1.516 

1.779 

7 370E48 

Wirmga UV14 + UVII 

1 04E15 

1.400 

1.400 

1.562 

O.OOOEOO 




1 420 

1 578 

3.050E48 




1.443 

1 598 

4.538E48 




1.467 

1.618 

5 655E48 




1 512 

1 656 

7.291E48 

Chiral Sigma Model 

4 06i;i4 

1 400 

1 400 

1.521 

0 OOOEOO 




1.434 

1.549 

3.973E48 




1 472 

1.581 

5 760£48 




1.505 

1.609 

6.976E48 




1 539 

1 638 

8.031E48 


outlined above. We assume that for a central density that gives a gravitational mass of 
14 Mq in the non-rotating case, / = Q = 0, and we then vary Q from zero to 12, The 
results are presented m Table 3; these should be compared with Table 1 of KW. 

From Table 3, we see that for each equation of state, a plot of Mg vs M’ gives an 
almost constant slope We solve equation (38) using the boundary condition J = 0 for 
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Figure 2. Evolution of the dimensionless angular momentum,; = J/M^, of a neutron star with 
the accretion of mass, for six representative equations of state. The dashed curves correspond to 
the approach of KW, and the solid curves to our more self-consistent approach 
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M' = 14Mq configuration We integrate it either until M' = where is 
mass corresponding to or until J = (the last entry for J in Table 3) whichc 
comes first 

The results of our computations for the angular momentum are illustrated in Fi| 
The x-axis corresponds to M'/M^ and the y-axis to J/M'^ The dashed cui 
correspond to the angular momentum evolution computed in the fashion of KW, w 
the solid curves represent the results of the more self-consistent (numerical) estim£ 
made by us as described above. The graphs indicate a slower rate of evolution of 
accreting neutron star's angular momentum in the rotational case as compared to 
case where rotational effects are not treated 

There are two factors contributing to the slower evolution of the angular moment 
that we obtain The less significant one is a comparatively smaller value of dAf^/d. 
The more significant factor is that the angular momentum T of the accreted partich 
systematically less in the rotational case. Our calculations include the fact that infall 
matenal co-rotatmg with the neutron star has a smaller cross-section for accretion th 
does counter-rotating infalling matter, so that the effective angular momentum 
reduced (eg, Misner, Thome, & Wheeler 1974) Within the realm of overlap of c 
approximations, the results we find agree reasonably well with those of Cook et 
(1994a, b) for accretion induced changes of angular momentum 


5. Discussion 

In presenting the results in the previous section, we have made the implicit assumptu 
that the magnetic field of the neutron star is too small to affect accretion Clear 
a quantitative estimate of this limit is in order The Alfven radius (r^), is defined by tl 
relationship (see Lamb, Pethick & Pines 1973) 

^ = (3 

where p and v are respectively the density and radial velocity in the accretion dis 
determines the location at which magnetic pressure channels the flow from a disk int 
an accretion column structure above the magnetic poles. Lamb, Pethick & Pines ( 1 97 
show that 




<2-6x 10® 




cm. 


(4( 


where = L 37 is the total luminosity in units of 10 ^^ ergs" 

= i^/ 10 ®cm and Bq is the magnetic field on the surface of the neutron star in gaus! 
The condition that r^<R implies that (taking M = 1-4 Mq and R^ = 1) 

B^<55 X 10'^ (41 


and is necessary for the scenario we have discussed to be fully self-consistent, but field 
somewhat higher than this value will not greatly modify our conclusions 
In our notation, L = {Fjy + Es)Mc^, with M the accretion rate According to ou 
calculations, typical values for + E'^) are of the order of 0 2 The luminosity L 37 = 



Disk Luminosity and Angular Momentum 


373 


would then correspond to an accretion rate -^5-6 x Such accretion rates are 

close to the ones estimated in X-ray binaries (Ghosh & Lamb 1991), so that our 
computations are relevant for systems with significant accretion on to old neutron stars 
whose surface magnetic fields have undergone substantial decay (to about 10® C). 

Under these circumstances of weak neutron star magnetic fields, we have shown that 
an incorporation of rotational effects always increases the disk luminosity, usually 
decreases the boundary layer luminosity, and always reduces the rate at which the 
neutron star’s angular momentum rises with accreted mass These effects are large 
enough to ment their consideration in analyses of observations of low-mass X-ray 
binaries. 
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Abstract. An accurate numencal survey of the density profiles cor- 
responding to the crusts of neutron stars for representative equation of state 
models is presented This will find application in calculations of thermal 
and magnetic evolution of neutron stars 
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1. Introduction 

The crust of a neutron star refers to a region, near the surface, that is composed of 
mainly neutron-nch bare nuclei which interact electrostatically and are arranged in 
a lattice The outer crust has small admixtures of free electrons and the inner crust has, 
in addition to electrons, ‘drip’ neutrons populating continuum states and permeating 
the nuclear lattice space. With increasing depth from the surface, and hence with 
increasing density, the neutron-to-proton ratio of the nuclei increases till a point is 
reached where the nuclei ‘dissolve’ and the composition is largely a homogeneous sea of 
neutrons This point is charactenzed as the crust bottom Although the crusts compnse 
a relatively small fraction of the total radius of the neutron star (usually about 10%), 
their importance derives from the fact that these possess interesting transport proper- 
ties, and are believed to play a central role in a variety of directly and indirectly 
observable properties of neutron stars. The former include glitches, X-ray emission, 
and the latter concern the evolution of neutron star magnetic fields (and hence the 
evolution of pulsars) and possible changes m the neutron star structure and composi- 
tion due to long-term accretion from a binary companion star In order to understand 
all these and related phenomena quantitatively, it is imperative to have accurate 
knowledge of the radial distribution of the density in the crust 
Calculations suggest that neutron stars possess a more or less flat density profile 
from the centre to the crust bottom (Arnett & Bowers 1977; Datta 1988), and that the 
crust contains negligible mass in comparison to the mass in the rest of the star The 
crust bottom (which is the boundary layer between the inner crust and the interior 
region) according to realistic calculations, is roughly equal to the equilibnum nuclear 
matter density, 2.4 x lO^^^gcm"^. This is about an order of magnitude less than the 
central densities obtained in neutron stars Across the crust, however, the density drop 
Ls extremely sharp - about fourteen orders of magnitude. Clearly, an accurate knowl- 
edge of the crustal density profile needs an accurate numerical treatment while 
integrating the equilibrium structure equations in this region of extremely rapid density 
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variation The equation of state (EOS) of neutron star matter is the key input part 
in all these calculations, the crustal properties depending sensitively on which pj 
lar EOS model is adopted to describe the neutron star interior In general, the stif 
EOS, the more the extent of the crust. In this paper, we present a numerical sur 
crustal density profiles for several neutron star configurations, and for a represen 
choice of the EOS models 


2. Equations of state 

The structure of neutron star depends sensitively on the EOS at high den; 
Although the main composition of degenerate matter at densities > 2.4 x 10^‘^g c 
that characterize neutron star interiors, is expected to be dominated by neut 
significant admixtures of other elementary particles (such as pions, muons 
hyperons) are not ruled out. A persistent problem in determining the EOS (nai 
pressure as a function of density) of such high density matter is what to choose fo 
interaction among the various particles, for which reliable experimental informati 
not available All calculations involve either extrapolations from known nuclear m- 
properties or field theoretical approaches Another unresolved problem is what m 
body technique is adequate for the purpose of evaluating the higher order correla 
terms for the pressure In this paper we do not address these problems, but choos< 
our purpose, eight EOS models based on representative neutron star matter mterac 
models A bnef description of these is given below 

A. Pandharipande (neutron matter): Pandhanpande (1971a) studied behavioui 
dense neutron matter usmg a many-body theory based upon the variational appro 
suggested by Jastrow (1955). The two-body wave function was taken as satisfy 
a simplified form of the Bethe-Goldstone equation, m which, terms representing 
Pauli exclusion principle were omitted but simulated by imposing a ‘healing’ constra 
on the wave function Instead of a state-dependent correlation function, an average \ 
used with spin dependence. To describe the nuclear interaction, the Reid nuclei 
nucleon potential was used (Reid 1968). This potential model is now considered to 
inadequate as it does not fully describe known nuclear matter properties 

R Pandbaripande (byperonic matter): At high densities that prevail m neutron si 
mtenors, it is possible that elementary particles heavier than neutrons, such as hypero] 
may also be present. The suggestion that hyperons may be additional baryonic constil 
ents of neutron star mtenors was first made by Ambartsumyan & Saakyan (1960) bas- 
on energetic arguments. Although our knowledge of hyperomc interactions is meag 
(pnmanly because of lack of expenmental data), there have been several theoretic 
attempts aimed at denvmg the EOS of baryonic hquid made up of neutrons, protor 
hyperons (A,Lo) One of the early such attempts is due to Pandhanpande (1971b), wt 
assumed the hyperomc potentials to be similar to the nucleon-nucleon potentials, bi 
altered suitably to represent the different isospin states The many-body method adopte 
to denve the EOS was similar as m the case A. Several EOS models of hyperomc matt^ 
have since been proposed by different authors. However, this still remains an ope 
problem because our knowledge of hyperon-nucleon and hyperon-hyperon mterac 
tions and their coupling constants have large uncertainties 
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C Bethe-Jobnson: Bethe & Johnson (1974) devised phenomenological potentials for 
nucleon-nucleon interaction that have more realistic short-range behaviour than the 
Reid potentials These authors then used the lowest order constrained variational 
method as given by Pandhanpande (1971 a) to calculate the EOS of neutron star matter. 
The work of Bethe & Johnson (1974) consists of two different parts: (1) detenmnation of 
the EOS for a pure neutron gas and (2) denvation of a hyperonic equation of state For 
purpose of illustration here, we have chosen their EOS model V (neutron matter), 

D. Waleckai The exchange of vector and scalar mesons among nucleons is known to 
provide the short-range repulsion and intermediate-range attraction in the nucleon- 
nucleon potential. The effective interaction will be characterized by the meson para- 
meters, such as their masses and coupling constants. If these parameters are specified, 
the meson wave functions can be easily solved as a function of density (assuming the 
wave functions to be space- time independent). From this, one can evaluate the stress 
tensor, which then provides the EOS. Walecka (1974) chose the meson parameters so as 
to reproduce the binding energy of nuclear matter (assumed to be — 15.75 MeV/nucleon) 
at a saturation density of 0 193 nucleons fm " More recent analysis of experimental data 
suggests the nuclear matter binding energy to be — 16.3 MeV/nucleon and a saturation 
density of 0153 nucleons fm"^ (Moller etal 1988) The EOS model of Walecka 
corresponds to pure neutron matter Both the scalar and vector mesons are assumed to 
have zero isospin by Walecka, so the model does not give the right asymmetry energy 
for nuclear matter. 

E. Friedmao-Pandharipande (neutrons): This model (Friedman & Pandharipande 
1981), is also based on the lowest order constrained vanational method. However, 
instead of the Reid potentials, the authors have made use of an improved pheno- 
menological nucleon-nucleon potential containing a two-body part as well as three- 
body correlations This interaction fits the nucleon-nucleon scattering cross-section 
data, deuteron properties, and also the nuclear matter properties rather well. 

F. Wiringa, Fiks & Fabrocini (1988): These authors gave a model of EOS for dense 
nuclear and neutron matter which includes three-nucleon interactions. This is a non- 
relativistic approach based on the variational method. The three-body potential 
considered by the authors includes long-range repulsive parts that are adjusted to give 
light nuclei binding energies and nuclear matter saturation properties. This work 
represents an improvement over the calculation of Fnedman & Pandharipande (1981) 
regarding the long-range attraction term in the Hamiltonian. The authors have given 
three models We consider here their model for beta-stable case* UV14 + UVII 
(neutrons, protons, electrons and muons). 

G. Prakasbf Ainsworth & Lattuner (1988) proposed EOS models for neutron stars, 
based on an extrapolation of the energy per particle of symmetric nuclear matter. The 
approach is empincal, and involves several parameters. The total energies are obtained 
by inclusion of the nucleon kinetic energy and the effect of finite forces between nucleons 
Four models have been suggested by the authors, we use model 1 for our purpose. 

K Sahu, Basu &Datta (1993) gave a field theoretical EOS for neutron-nch matter in 
beta equihbnum based on the chiral sigma model. The model includes an isoscalar vector 
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field generated dynamically and reproduces the empirical values of the nuclear 
saturation density and bmding energy and also the isospin symmetry coefficie 
asymmetric nuclear matter. The energy per nucleon of nuclear matter according to 
Basu & Datta (1993) is in very good agreement, up to about four times the equili 
nuclear matter density, with estimates inferred from heavy-ion collision expenmenta 
Models A, B and G are relatively soft EOS, models E, C and F are roughly interm 
m theu: stif&iess whereas models D and H are very stiff EOS. The composite EOS f 
entire span of neutron star densities was constructed by joining the selected high d 
EOS to that of Negele & Vauthenn (1973) for the density range (10^"^ — 5 x 10^°)g( 
Baym, Pethick & Sutherland ( 1 97 1 ) for densities down to 1 0^ g cm “ ^ and Feyr 

Metropolis & Teller (1949) for densities less than 10^gcm“^. 


3. Calculations 

The Tolman-Oppenheimer-Volkoff (Oppenheimer & Volkoff 1939) equations that 
be integrated numerically to obtain stable, non-rotating neutron star structure a 

dr 

dP _ G{m + 4Tcr^Pjc^)(p + P/c^) 
r{r-2Gm/c^) ’ 



Figure 1. Pressure as a function of density for the EOS models (A-H) 
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Table 1. Neutron star structure parameters for the EOS Models (A-H) 


EOS 

M/Mq 

R (km) 

Pc (gem'®) 

A (km) 

Mcr/M 

Pandbaripande 

1,330 

10 006 

1 669i;+ 15 

0 963 

3 024£ - 02 

(neutrons) 

1 400 

09 867 

1 S30E + 15 

0 858 

2 537 E - 02 


1.658 

08 492 

4 100£; + 15 

0.404 

8 03l£-03 

Pandhanpande 

1 330 

07 921 

3 606£+15 

0 484 

9 189£ - 03 

(hyperons) 

1 400 

07 437 

4 795£: + 15 

0 361 

5 863E - 03 


1.414 

07 082 

6 017£:+15 

0 300 

4 426£; - 03 

Bethe- Johnson 

1.330 

10 782 

1 283E + 15 

1 149 

3 945f: - 02 

Model V (neutrons) 

1.400 

10 691 

1.373£ + 15 

1 044 

3 41l£-02 


1.758 

09 200 

3.466£ + 15 

0 454 

9 565£ - 03 

Walecka 

1 330 

12.252 

7.158^;+ 14 

1 442 

5 482£ - 02 


1 400 

12 280 

7 415i; + 14 

1 387 

5 415f;-Q2 


2 285 

11 216 

2 1745 + 15 

0 452 

1.040£ - 02 

Fnedman- Pandbaripande 

1.330 

10 824 

1 332E+15 

1 134 

3 859E-02 

(neutrons) 

1 400 

10.701 

1 413E+15 

1 023 

3.297E-02 


1 986 

09 028 

3.450E+15 

0.307 

5 531E-03 

Wiringa et al 

1 330 

11 122 

1.000E+ 15 

1.120 

3.543f; - 02 

(UV14 -b UVII) 

1 400 

11 115 

1 043£;+15 

1.039 

3.167J? - 02 


2 189 

09.800 

2.830f?+lS 

0.294 

4.946i; - 03 

Prakash, Ainsworth 

1 330 

11 384 

1.507E+ 15 

1.702 

7 376£: - 02 

L Lattimer Model 1 

1,400 

10 987 

1.808E+ 15 

1.456 

5 781E: - 02 


1 510 

09 483 

3 590i; + 15 

0.901 

2.744£^ - 02 

Sahu, Basu 

1 330 

14 685 

3 952E + 14 

3.128 

2 518£ - 01 

Datta 

1.400 

14 761 

4.063E + 14 

2.970 

2 342f; - 01 


2 593 

14 096 

1 345£ + 15 

0.975 

4 588f: - 02 


NOTE* In Tables 1-9, the numbers following the letter E represent powers of ten. 


where P and p stand for pressure and density and m is the gravitational mass contained 
within a sphere of radius r. The integration of equations (1) and (2) is done for a given 
EOS which specifies P(p\ and with initial condition p = Pc, a chosen value for the 
central density. The numerical integration started from the centre of the star is 
terminated at the radius of the star, defined as the point where the density 

7 86gcm“^, expected to correspond to the neutron star surface The total gravita- 
tional mass is then given by M = m(R). 

We integrate equations (1) and (2) simultaneously using Runge-Kutta fourth order 
method, having an adaptive step size control The step size is adjusted for the required 
accuracy from a calculation of the local truncation error (Antia 1991) We find that this 
provides an extremely accurate method to integrate equations (1) and (2) for neutron 
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Table 2. Crustal density profiles Pandharipande (neutrons) 


M = 1 330 Mq 

M = 1 400 Mq 

M = 1.658 M 

R = 10 006 (km) 

R = 9.867 (km) 

R = 8 492 (kj 

(1 - r/R) 

(gem 

(1 - r/R) 


(1 - r/R) 

P 

(g Cl 

O.OOOOOi; + 00 

7.8665 + 00 

0 000005 + 00 

7.8635 + 00 

0.000005 + 00 

7 866i 

1 13773f; - 09 

1.3195+01 

1 288915 - 09 

1.3605 + 01 

9.917955 - 10 

1 447i 

2.39423£ - 08 

4 9725 + 01 

2 230905 - 08 

51205 + 01 

1.275795 - 08 

5 428i 

1.41172jE: - 07 

1 8745+02 

1.310955 - 07 

1 9275 + 02 

7 444055 - 08 

2.037f 

6 75S15f; - 07 

7 0665 + 02 

6.218655 - 07 

7 2525 + 02 

3.468165 - 07 

7.642£ 

2 35836i; - 06 

2 6645 + 03 

2.170805 - 06 

2 7305 + 03 

1.197435 - 06 

2.8675 

6 95373i: - 06 

1.0045 + 04 

6.340945 - 06 

1.0275 + 04 

3 47397i; - 06 

1 0765 

1.94642i; - 05 

3.7865 + 04 

1.770385 - 05 

3.8665 + 04 

9.739525 - 06 

4 0375 

5 USllE - 05 

1 4275 + 05 

4 671515 - 05 

1 4555 + 05 

2.537655 - 05 

1 5155 

1.26075i; - 04 

5.3825 + 05 

1.141605 - 04 

5.4775 + 05 

6.238155 - 05 

5 6835 

2 88548f; - 04 

2 0295 + 06 

2 608805 - 04 

2 0615 + 06 

1.404455 - 04 

2 1325 

6 06912£’ - 04 

7,6495 + 06 

5.471535 - 04 

7 7585 + 06 

2.932735 - 04 

8 0015 

1 14403f; - 03 

2.8845 + 07 

1.030125 - 03 

2 9205 + 07 

5.562135 - 04 

3 0025 

2.05310£ - 03 

1.0875 + 08 

1.847385 - 03 

1 0995 + 08 

9 865455 - 04 

1.1265 

3 46132£: - 03 

4.0995 + 08 

3.112335 - 03 

4.1365 + 08 

1 659945 - 03 

4.2265 

5 441 70£; - 03 

1 5455 + 09 

4 891285 - 03 

1-.5575 + 09 

2.610385 - 03 

1.5865 

8.57267£ - 03 

5 8265 + 09 

7 700805 - 03 

5 8595 + 09 

4.107925 - 03 

5 9505 - 

1 28894£ - 02 

2 1965 + 10 

1.156915 - 02 

2.2055 + 10 

6.207215 - 03 

2.2335 - 

1.89605f: - 02 

8.2805 + 10 

1.702605 - 02 

8.2995+10 

9 140445 - 03 

8.3775 - 

2 54733£ - 02 

3 1225 + 11 

2.287865 - 02 

3.1245+ 11 

1.228325 - 02 

3.1435 - 

3 04385£ - 02 

1 1775 + 12 

2.734825 - 02 

1.1765 + 12 

1 470045 - 02 

1.1795 H 

3 394945 - 02 

4 4375 + 12 

3 051035 - 02 

4.4255+12 

1.641645 - 02 

4 4255 H 

4 029995 - 02 

1 6735 + 13 

3 622395 - 02 

1.6655 + 13 

1.946755 - 02 

1 6605 ^ 

5 348545 - 02 

6 3065 + 13 

4 810265 - 02 

6 2685 + 13 

2.597685 - 02 

6 2305 -i 

9 626145 - 02 

2 4005 + 14 

8.696105 - 02 

2.4005 + 14 

4.753095 - 02 

2 4005 4 


Note. In the crustal density profiles presented here and in the subsequent tables, the th 
mass column corresponds to that of the maximum stable mass. 


Star configurations varying widely in their density profiles. The method, of course, giv 
bulk properties such as masses, radii and moments of mertia that agree with the valu 
already reported in the literature. 
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Table 3. Crustal density profiles Pandharipande (hyperons) 


M = 1330 Mo 

M = 1 400 Mo 

M = 1 658 Mo 

R = 7 921 (km) 

R = 7 437 (km) 

R = 7 082 (km) 

(1 - r/R) 

(g cm-®) 

(1 - r/R) 

(gem 

(1 - r/R) 

(g cm ®) 

0 00000£ + 00 

7 8691? + 00 

0 000001? + 00 

7 8731? + 00 

0 000001? + 00 

7.8721 + 00 

1 07340f: - 09 

1 3861? + 01 

5 541091? - 10 

1.3101: + 01 

7.304791? - 10 

1.3861 + 01 

1 67489i; - 08 

5 2121? + 01 

1 222701? - 08 

4 9371? + 01 

1 143051' - 08 

5.2061 + 01 

9.81838.B - 08 

1 9591' + 02 

7 205071? - 08 

1 8611? + 02 

6 693611 - 08 

1.9561 + 02 

4 63059i; - 07 

7 3661? + 02 

3 448561? - 07 

7 0131? + 02 

3 153451 - 07 

7.3491 + 02 

1 61124i; - 06 

2 7691? + 03 

1.215811? -06 

2.6431? + 03 

1 098901 - 06 

2 7611 + 03 

4 69519E - 06 

1.0411? + 04 

3 557921? - 06 

9.9611? + 03 

3.197851 - 06 

1.0371 + 04 

1 3188l£' - 05 

3 9141? + 04 

1 003631? - 05 

3.7541? + 04 

8 989331 - 06 

3 8981 + 04 

3 44931i; - 05 

1 4711? + 05 

2 632821? - 05 

1.4151? + 05 

2 347701 - 05 

1.4641 + 05 

8 494831? - 05 

5 5311? + 05 

6 512091? - 05 

5.3331? + 05 

5 787401 - 05 

5.5021 + 05 

1 920581? - 04 

2 0791? + 06 

1 476001? - 04 

2.0101? + 06 

1 305951 - 04 

2 0671 + 06 

4.024301? - 04 

7 8171? + 06 

3.114791? - 04 

7.5751? + 06 

2.742081 - 04 

7 7671 + 06 

7 633181? - 04 

2.939.B + 07 

5 915711? - 04 

2.8551? + 07 

5.203081 - 04 

2.9181 + 07 

1 356301? - 03 

1.1051? + 08 

1.051991? - 03 

1.0761? + 08 

9.301931 - 04 

1 0961 + 08 

2 285911? - 03 

4.153£ + 08 

1.778141? - 03 

4.0561? + 08 

1.557241 - 03 

4.1201 + 08 

3 593961? - 03 

1.5611? + 09 

2.800541? - 03 

1.5291? + 09 

2.450391 - 03 

1 5481 + 09 

5 660661? - 03 

5 8701? + 09 

4.423191? - 03 

5 7611? + 09 

3.865071 - 03 

5.8151 + 09 

8 555491? - 03 

2.2071? + 10 

6.701781? - 03 

2.1711? + 10 

5.850491 - 03 

2 1851+10 

1.258521? - 02 

8 2961? + 10 

9.870521? - 03 

8.1831? + 10 

8.614291 - 03 

8 2091 + 10 

1 708221? - 02 

3.1191? + 11 

1.344601? -02 

3 0841? + 11 

1.172951 -02 

3.0841 + 11 

2.033231? - 02 

1.1731? + 12 

1 602611? - 02 

1 1621? + 12 

1.398061 - 02 

1 1591 + 12 

2 268431? - 02 

4 4081? + 12 

1 789471? - 02 

4 3811? + 12 

1.561061 - 02 

4.3541 + 12 

2 687351? - 02 

1 6571? + 13 

2 120471? - 02 

1.6511? + 13 

1 847891 - 02 

1.6361 + 13 

3 57424.E - 02 

6.2301? + 13 

2 827031? - 02 

6.2231? + 13 

2 462111 - 02 

6.1461 + 13 

6 111921? -02 

2 4001? + 14 

4.849711? - 02 

2 4001? + 14 

4.241781 - 02 

2.4001 + 14 


4. Results and discussions 

rhe EOS models (A-H) are illustrated in Fig. 1 (which is a plot of pressure versus 
iensity) This plot refers to the high density region corresponding to the neutron star 
ntenor. We have already mentioned in section 2 the EOS that is used for the crustal 
egion. It is important to note that because of substantial variation in the various EOS 
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models (A-H), the crustal thicknesses and density profiles are expected to exh 
vanation even though the crustal EOS is assumed to be the same for all cases. 

In Table 1, we have summarized the results for the mass (M), radius (R\ cen 
density (pj, the crustal thickness (A) and the mass contained m the crust as a fractioi 
the total gravitational mass We have done this for three configurations of 


Table 4. Crustal density profiles: Bethe- Johnson Model V (neutrons). 


M = 1.330 Mq 

M = 1 400 Mq 

M = 1.658 Mq 

R = 10.782 (km) 

R = 10.691 (km) 

R = 9 200 (km) 

(1 - rlR) 

(g cm“®) 

(1-r/R) 

(g cm"^) 

(1 - r/R) 

P 

(g cm"^ 

0 OOOOOf? + 00 

7.8655 + 00 

0 000005 + 00 

7 8655 + 00 

0 000005 + 00 

7 8705 + 

9 62517f: - 10 

1 2745 + 01 

5 121155 - 10 

1 1935 + 01 

1.181525 - 09 

1 4835 + 

2.584 71£ - 08 

4 8115 + 01 

2.158505 - 08 

4.5175 + 01 

1 385935 - 08 

5.5615 + 

1.52912i;-07 

1 8175 + 02 

1.284315 - 07 

1 7115 + 02 

8 067565 - 08 

2.0845 + 

l.mSiE - 07 

6 8595 + 02 

6.315125-07 

6 4795 + 02 

3.734275 - 07 

7.8145 + 

2 60124i; - 06 

2.5905 + 03 

2.254725 - 06 

2 4545 + 03 

1.282365 - 06 

2 9295 + 1 

7.68367i; - 06 

9 7805 + 03 

6.704875 - 06 

9.2925 + 03 

3.713895 - 06 

1 0985 + ( 

2 15801£ - 05 

3 6935 + 04 

1 889685 - 05 

3 5195 + 04 

1.039455 - 05 

4 1165 + C 

5 71082£; - 05 

1 3945 + 05 

5 025365 - 05 

1 3335 + 05 

2.704605 - 05 

1.5435 + C 

1.404775 - 04 

5.2655 + 05 

1.240645 - 04 

5 0475 + 05 

6.637165 - 05 

5.7845 + C 

3 217955 - 04 

1 9885 + 06 

2.855965 - 04 

1.9115 + 06 

1 492435 - 04 

2 1685 + 0 

6.789335 - 04 

7.5075 + 06 

6 069575 - 04 

7.2385 + 06 

3 106765 - 04 

8.1295 + 0 

1.282515 - 03 

2 8355 + 07 

1 147195 - 03 

2 7415 + 07 

5 889685 - 04 

3 0475 + 0 

2 302145 - 03 

1.0705 + 08 

2.065315 - 03 

1.0385 + 08 

1 042515 - 03 

1.1425 + 0 

3.883115 - 03 

4.0415 + 08 

3.492725 - 03 

3.9315 + 08 

1 754885 - 03 

4.2825 + 0. 

6 105975 - 03 

1.5265 + 09 

5.502945 - 03 

1 4895 + 09 

2.758585 - 03 

1.6055 + 0) 

9.576205 - 03 

5.7625 + 09 

8.643905 - 03 

5 6385 + 09 

4.334855 - 03 

6 0175 + 0) 

1 449655 - 02 

2 1765 + 10 

1.312035 - 02 

2 1355 + 10 

6 543685 - 03 

2.2565 + It 

2.129805 - 02 

8 2155 + 10 

1 929675 - 02 

8 0865 + 10 

9 636275 - 03 

8 4565 + 1C 

2.818675 - 02 

3.1025 + 11 

2.562955 - 02 

3 0625 + 11 

1.272525 - 02 

3.1705 + 11 

3.388685 - 02 

1.1715 + 12 

3.085825 - 02 

1.1605+12 

1 531705 - 02 

1 1885 + 12 

3 782175 - 02 

4.4235 + 12 

3.446675 - 02 

4.3925 + 12 

1.711855 - 02 

4.4545 + 12 

4.491995 - 02 

1.6705 + 13 

4 096005 - 02 

1 6635 + 13 

2.033585 - 02 

1.6705 + 13 

5 927135 - 02 

6.3065 + 13 

5.410595 - 02 

6.2995 + 13 

2 718255 - 02 

6.2595 + 13 

1 065395 - 01 

2.4005 + 14 

9.761665 - 02 

2.4005 + 14 

4.937175 - 02 

2.4005 + 14 
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Table 5. Crustal density profiles Walecka 


M = 1 330 Me 

M = 1 400 Mq 

M = 1 658 Mo 

R = 12.252 (km) 

iZ = 12 280 (km) 

R = 11.216 (km) 

(1 - r/R) p 

(1 - rfR) p 

(1 - r/5) p 

(g cm ^) 

(g cm-^) 

(g 


0.0000015 + 00 

1 23191£' - 09 
3.16611i;-08 
1.86840E-07 
9.00051£;-07 

3 16390£: - 06 
9.32S44i; - 06 

2 61696£ - 05 

6 91322i; - 05 

1 70007£ - 04 

3.88916£ - 04 
8.20653£ - 04 
1.5496(r£; - 03 

2 77S95E - 03 
4.68399f; - 03 

7.30487jE; - 03 
1.15413£^ - 02 
1.74576f; - 02 
2.55827f; - 02 
3 . 44452 ^; - 02 

4.11286E-02 
4.58221i; - 02 
5.4115615-02 

7 09277jE: - 02 
1.17734i; - 01 


7 86115 + 00 
1 285E + 01 
4.842i: + 01 

1 825 £; + 02 
6 . 877 ^; + 02 

2 592£ + 03 
9 767 E + 03 
3.681£ + 04 
1.387f; + 05 
5.228£ + 05 

1 970E + 06 

7 426£ + 06 
2.799S + 07 
1 055jE + 08 
3.975i; + 08 

1.498£ + 09 
5.646£ + 09 
2.128i; + 10 

8 019£ + 10 

3 02215+ 11 

1 139i; + 12 

4 292f; + 12 
1 618£ + 13 
6 097E + 13 
2.400i; + 14 


0 OOOOOE + 00 
8.38888£; - 10 
2.80625£; - 08 
1.66453f; - 07 
8 10853£; - 07 

2 87466£; - 06 
8 51191£; - 06 
2.39565£' - 05 
6.3S046f; - 05 

1 56575f; - 04 

3.59261£; - 04 
7.60393f? - 04 
1.43749E - 03 

2 58236£ - 03 
4.35895f; - 03 

6.83400£; - 03 
1 07640f; - 02 

1 62969^; - 02 

2 39257E - 02 

3 21742f; - 02 

3 84200i; - 02 
4.28220f; - 02 
5.07298E - 02 
6.67522E - 02 
1.12973£-01 


7 863i; + 00 
1 236 ^; + 01 

4 672E + 01 

1 766£ + 02 

6 679E + 02 

2 525£! + 03 
9.546£ + 03 

3 609£; + 04 
1.3645 + 05 

5 1585 + 05 

1 9505 + 06 

7 3735 + 06 

2 7885 + 07 
1.0545 + 08 

3 9845 + 08 

1 5065 + 09 

5 6955 + 09 
2.1535+10 

8 1405 + 10 
3 0785 + 11 

1 1645 + 12 
4.3995+12 
1.6635 + 13 

6 2875 + 13 

2 4005 + 14 


0 000005 + 00 
8.180355 - 12 
5.395875 - 09 

3 348225 - 08 

1 817475 - 07 

7 108995-07 

2 184965 - 06 

6 342565 - 06 
1 716505 - 05 

4 384725 - 05 

1 034775 - 04 

2 270395 - 04 
4.364435 - 04 

7 978675 - 04 

1 364795 - 03 

2 199645 - 03 
3.464605 - 03 

5 323325 - 03 
7 894485 - 03 
1 090135 - 02 

1.314755 - 02 
1 473595 - 02 

1 744905 - 02 

2 336655 - 02 
4.034365 - 02 


7.8735 + 00 
8 1045 + 00 
3.1195 + 01 
1.2005 + 02 
4 6215 + 02 

1 7785 + 03 
6 8455 + 03 
2.6355 + 04 
1 0145 + 05 
3.9035 + 05 

1.5025 + 06 
5.7825 + 06 
2.2255 + 07 
8.5655 + 07 
3.2975 + 08 

1.2695 + 09 
4.8845 + 09 
1.8805 + 10 
7.2355 + 10 
2.7855 + 11 

1.0725 + 12 
4.1255 + 12 
1.5885 + 13 
6.1115 + 13 
2.4005 + 14 


gravitational mass, namely, 1.33 Mq, 1 4 Mq and the maximum stable mass. A look at 
Table 1 shows that, for the EOS models considered here, the maximum neutron star 
mass (M) and their corresponding radu (R) he in the range (1 414-2 593) Mq, and 
(7.082-14.096) km respectively, with the lower value of each of the quantities corre- 
sponding to the Pandhanpande hyperomc matter EOS (a very soft one) and the higher 
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to the EOS due to Sahu, Basu & Datta (a very stiff EOS). The values of the 
thickness (A) for maximum mass configurations lie in the range (0 294-0 975 
Clearly, the crust thickness, and therefore the crustal density profile, depend sensit 
on the EOS model chosen to describe neutron star matter Of the three i 
configurations listed in Table 1, the 1 4 Mq neutron star models are of partic 
interest because analysis of available binary pulsar data suggests the pulsar mass< 


Table 6. Crustal density profiles. Friedman- Pandhaxipande (neutrons) 


M = 1.330 Mo 

M = 1.400 Mo 

M = 1.658 Me 

JZ = 10 824 (km) 

R = 10 701 (km) 

R = 9.028 (km) 

(1 - r/R) 

(g 

(1 - r/R) 

(g cm-®) 

' (l-r/5) 

P 

(g cm- 

0.00000£ + 00 

7 860jE; + 00 

0 00000£ + 00 

7.8635 + 00 

0 000005 + 00 

7 8615 + 

8 52656i; ~ 10 

1 255E + 01 

1 37922f; - 09 

1 3505 + 01 

1 298135 - 09 

1 6735 + 

2.55103£'-08 

4.744£ + 01 

2.53386£' - 08 

5>0865 + 01 

1.129135 - 08 

6 2375 + 

l.S1130i;-07 

1.793£ + 02 

1 49079i; - 07 

1 9155 + 02 

6 476145 - 08 

2 3255 + 

7.33104f? - 07 

6 773E + 02 

7.08975£ - 07 

7.2135 + 02 

2.898985 - 07 

8 6675 + 

2.58950£: - 06 

2 559£ + 03 

2.46504i; - 06 

2 7165 + 03 

9.72901F; - 07 

3 23l5 + ( 

7.65971i; - 06 

9 671f; + 03 

7 24485^: - 06 

1.0235 + 04 

2 788385 - 06 

1 2045 + ( 

2 15325£' - OS 

3 654J5; + 04 

2.02704i; - 05 

3.8525 + 04 

7 749265 - 06 

4.4905 + C 

5 70408F; - 05 

1 381jE; + 05 

5 34425i; - 05 

1.4515 + 05 

1 997665 - 05 

1.6745 + 0 

1.40445f: - 04 

5.217i;+05 

1 30935f; - 04 

5 4635 + 05 

4.870085 - 05 

6.2395 + 0 

3 22040.E - 04 

1.971i; + 06 

2 98868f; - 04 

2 0575 + 06 

1 094015 - 04 

2.3265 + 0i 

6 80S00£ - 04 

7 449jE; + 06 

6 26729i; - 04 

7.7485 + 06 

2.233415 - 04 

8.6705 + 0( 

1 28579£ - 03 

2 S15E + 07 

1 18248i; - 03 

2.9185 + 07 

4.226605 - 04 

3.2325 + O'i 

2 30933i: - 03 

1.064£ + 08 

2.11855i;-03 

1 0995 + 08 

7.507975 - 04 

1 2055 + 08 

3 89716i; - 03 

4.019£' + 08 

3.56727J? - 03 

4 1385 + 08 

1 247885 - 03 

4 4915 + 08 

6 130091; - 03 

1 519i; + 09 

5.60437.E - 03 

1.5585 + 09 

1 959885 - 03 

1.6745 + 09 

9.61827f; - 03 

5 7ZSE + 09 

8 77043i; - 03 

5.8695 + 09 

3 065895 - 03 

6 2415 + 09 

1 45658f; - 02 

2.168£: + 10 

1.326205 - 02 

2.2105 + 10 

4.617475 - 03 

2.3265 + 10 

2.14052£^ - 02 

8 193i; + 10 

1.949555 - 02 

8.3245 + 10 

6 803115-03 

8 6735 + 10 

2 87813JE: - 02 

3 096£: + 11 

2 615645 - 02 

3.1355 + 11 

9 091915 - 03 

3.2335 + 11 

3.43830£ - 02 

1 170£ + 12 

3.123895 - 02 

1.1805 + 12 

1.086895 - 02 

1.2055 + 12 

3.83380£: - 02 

4 420£; + 12 

3.483295 - 02 

4.4465 + 12 

1 213445 - 02 

4 4925 + 12 

4.54710£' - 02 

1.670f; + 13 

4 133805 - 02 

1 6745 + 13 

1.438775 - 02 

1 6755 + 13 

5 99013i; - 02 

6 311i; + 13 

5 445975 - 02 

6 3055 + 13 

1.921925 - 02 

6 2435+13 

1.04735^:- 01 

2.400i; + 14 

9 556925 - 02 

2 4005 + 14 

3.401625 - 02 

2 4005 + 14 
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Cable 7. Ciustal density ptoiiles Wiringa et al (UV14 + UVII). 


M = 1.330 Afo 

M = 1 400 Mq 

M = 1.658 Mq 

R = 11.122 (km) 

R = 11.115 (km) 

R = 9 800 (km) 

(1 - r/R) 

/> 

(g cm"®) 

(1 - rfR) 

1 

P \ 

(g cm"®) 

(1 - rfR) 

(gem ®) 

i.ooooof; + 00 

7 8655 + 00 

0 000005 + 00 

7.8635 + 00 

0 000005 + 00 

7 8625 + 00 

1 79318i: - 10 

9 5455 H- 00 

3.122465 - 11 

8.1235 + 00 

1 339835 - 09 

1 6845 + 01 

1 84269i; - 08 

3 6495 + 01 

1.351785 - 08 

3.1275 + 01 

1 105685 - 08 

6.2745 + 01 

1 HSOOi; - 07 

1 3955 + 02 

8.718345 - 08 

1.2045 + 02 

6 259455 - 08 

2 3385 + 02 

> 880075 - 07 

5.3325 + 02 

4.629225 - 07 

4.6365 + 02 

2.782945 - 07 

8 7105 + 02 

! 215765 - 06 

2.0385 + 03 

1.811265 - 06 

1 7855 + 03 

9 352525 - 07 

3 2465 + 03 

p 740565 - 06 

7.7925 + 03 

5 621995 - 06 

6 8725 + 03 

2.676725 - 06 

1 2095 + 04 

922715 - 05 

2.9795 + 04 

1.616285 - 05 

2 6465 + 04 

7.435465 - 06 

4.5065 + 04 

.190155 - 05 

1.1395 + 05 

4.415695 - 05 

1 0195 + 05 

1.915345 - 05 

1 6795 + 05 

.299665 - 04 

4 3535 + 05 

1.116435 - 04 

3 9225 + 05 

4 668045 - 05 

6.2565 + 05 

.037715 - 04 

1 6645 + 06 

2.641225 - 04 

1 5105 + 06 

1 048725 - 04 

2.3315 + 06 

.603095 - 04 

6.3615 + 06 

5.822895 - 04 

5 8135 + 06 

2.138575 - 04 

8.6865 + 06 

.253845 - 03 

2 4325 + 07 

1.111035-03 

2 2385 + 07 

4 046725 - 04 

3.2365 + 07 

277295 - 03 

9 2955 + 07 

2.033115-03 

8.6175 + 07 

7 187235 - 04 

1 2065 + 08 

.884605 - 03 

3 5535 + 08 

3.495935 - 03 

3.3185 + 08 

1 201735 - 03 

4 4935 + 08 

178315 - 03 

1 3585 + 09 

5.613925 - 03 

1.2775 + 09 

1 875295 - 03 

1.6745 + 09 

702515 - 03 

5.1935 + 09 

8.776125 - 03 

4.9185 + 09 

2.933505 - 03 

6 2385 + 09 

481275 - 02 

1 9855 + 10 

1.346255 - 02 

1.8935 + 10 

4.418145 - 03 

2 3245 + 10 

181865 - 02 

7.5885 + 10 

1.990085 - 02 

7.2905 + 10 

6.509135 - 03 

8 6615 + 10 

973005 - 02 

2.9015 + 11 

2 733595 - 02 

2 8075 + 11 

8 701285 - 03 

3.2275 + 11 

564095 - 02 

1.1095 + 12 

3.285365 - 02 

1.0815 + 12 

1 040315 - 02 

1.2035 + 12 

979735 - 02 

4 2395 + 12 

3 673445 - 02 

4 1605 + 12 

1 161475 - 02 

4 4815 + 12 

712005 - 02 

1 6205 + 13 

4.350515 - 02 

1.6025 + 13 

1.376625 - 02 

1 6705 + 13 

180975 - 02 

6 1955 + 13 

5 715565 - 02 

6 1675 + 13 

1.830395 - 02 

6.2215 + 13 

007415 - 01 

2 4005 + 14 

9.349535 - 02 

2.4005 + 14 

2 997435 - 02 

2.4005 + 14 


:se binanes to be close to this value For the 1 4Mo models, the A is in the range 
361-2 970) km, a much wider range than for the maximum mass cases. In Tables 2-9, 
give the crustal density variation with the depth m normalized units with respect to 
radius (1 — rjR), (where r is the radial distance from the centre of the star) from the 
face of the neutron star for the three mass configurations. These tables will be of 
ect application in calculations of magnetic field evolution of pulsars. The same thing, 
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namely, the crustal density variations, is graphically illustrated in Figs 2 and 
1 4 Mq configurations for the various EOS models The dashed lines in these fig 
represent the core-inner crust boundary (which corresponds to r = R — A) anc 
dot-cum-dashed lines represent the surface of the neutron star (r = R) The vanatii 
sharpest for the softest EOS considered here (Pandhanpande hyperomc matter) 
gets to be comparatively less sharp as the EOS gets to be more stiff. The dependence 


Table 8. Crustal density profiles: PraJcash, Ainsworth Lattimer Model 1 


M = 1 330 Mq 

M = 1.400 Mq 

M = 1.658 Mq 

R = 11.384 (km) 

R = 10 987 (km) 

= 9 483 (km; 

(1 - r/R) 

(gem 

(l-r/R) 

(gem 

(l-r/R) 

P 

(g cm’ 

0.00000£ + 00 

7 8645 + 00 

0 000005 + 00 

7.8615 + 00 

0 000005 + 00 

7.8685 

1.23204f; - 09 

1 2995 + 01 

5.643845 - 10 

1.2025 + 01 

7 204665 - 10 

1 2975 

2.87558f; - 08 

4 8935 + 01 

2 269685 - 08 

4.5395 + 01 

1 706005 - 08 

4 8935 - 

1 695iof? - or 

1 8435 + 02 

1.345505 - 07 

1 7145 + 02 

1.007395 - 07 

1 8465 - 

8 13995£ - 07 

6 9455 + 02 

6 592615 - 07 

6.4765 + 02 

4.840155 - 07 

6 9665 - 

2 85212^; - 06 

2 6165 + 03 

2.347115 - 06 

2 4465 + 03 

1.708215 - 06 

2 6285 - 

8.407995 - 06 

9.8575 + 03 

6.972315 - 06 

9.2385 + 03 

5.009585 - 06 

9 9175 - 

2.355635 - 05 

3.7135 + 04 

1.959995 - 05 

3.4895 + 04 

1 413275 - 05 

Z.74:2E - 

6.224945 - 05 

1.3995 + 05 

5.208635 - 05 

1.3185 + 05 

3.714745 - 05 

1.4125 i 

1 527715 - 04 

5 2705 + 05 

1 283045 - 04 

4.9785 + 05 

9.187525 - 05 

5 3275 H 

3 496395 - 04 

1 9865 + 06 

2.953965 - 04 

1.8805 + 06 

2.087105 - 04 

2.0105 H 

7.372035 - 04 

7.4805 + 06 

6 286985 - 04 

7.1025 + 06 

4.401245 - 04 

7 5855 4 

1 390125 - 03 

2 8185 + 07 

1 185485 - 03 

2 6835 + 07 

8.358405 - 04 

2 8625 4 

2 494185 - 03 

1 0625 + 08 

2.134885 - 03 

1.0135 + 08 

1.489475 - 03 

1 0805 4 

4 204875 - 03 

4 0005 + 08 

3.611335 - 03 

3 8275 + 08 

2.515145 - 03 

4 0745 4 

6 607105 - 03 

1 5075 + 09 

5 692285 - 03 

1 4465 + 09 

3 959055 - 03 

1.5375 + 

1.040935 - 02 

5 6775 + 09 

8 982415 - 03 

5 4605 + 09 

6.246935 - 03 

5 8015 + 

1.568005 - 02 

2 1395 + 10 

1 356855 - 02 

2 0625 + 10 

9 430945 - 03 

2 1895 + 

2.296565 - 02 

8.0575 + 10 

1 990825 - 02 

7.7905 + 10 

1.386315 - 02 

8.2585 + 

3 278495 - 02 

3.0365 + 11 

2 853765 - 02 

2.9425 + 11 

1.983625 - 02 

3.1165 + 

3.719075 - 02 

1.1445 + 12 

3.261605 - 02 

1.1115+12 

2.249555 - 02 

1 1765 + 

4 155135 - 02 

4.3085 + 12 

3.636845 - 02 

4.1985 + 12 

2 514235 - 02 

4.4365 + 

4 961695 - 02 

1 6235 + 13 

4 341545 - 02 

1.5865 + 13 

3.019745 - 02 

1 6745 + 

6 636625 - 02 

6 1155 + 13 

5 796965 - 02 

5 9895 + 13 

4.118305 - 02 

6.3165 + 

1.495365 - 01 

2.4005 + 14 

1.325275 - 01 

2.4005 + 14 

9.501825 - 02 

2.4005 + 
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Table 9. Crustal density profiles. Sahu, Basu & Datta 


M = 1 330 Me 

M = 1 400 Me 

M = 1 658 Mo 

R = 14.685 (km) 

R = 14.761 (km) 

R = 14.096 (km) 

(1 - r/R) p 

(1 - r/R) p 

(1 - r/R) p 

(g cm"®) 

(g cm-®) 

(g cm“®) 


0 00000 ^; + 00 

1 68177f; - 10 
2.51686£ - 08 
1.542025 - 07 

8 19531i; - 07 

3.13725£ - 06 
9.61781£; - 06 

2 75907f; - 05 
7 480S5f; - 05 

I. 88519£ - 04 

4 4246915 - 04 

9 61559f? - 04 
1 8438315 - 03 

3 35559£; - 03 

5 73697 E - 03 

9 15009£? - 03 
1 43445f; - 02 
2.1874315-02 
J 21388f? - 02 
t.41761i;-02 

5 25674f: - 02 
5 85730f? - 02 
}.90648i; - 02 
) 05533i? - 02 

J. 13000i; - 01 


7 861£; + 00 

8 839£ + 00 
3 391S + 01 
1 301£ + 02 
4.990i; + 02 

1 914£ + 03 

7 344f; + 03 

2 817.E + 04 
1 081f: + 05 
4.146f; + 05 

1 591f; + 06 
6.102f; + 06 

2 341i; + 07 

8 979E + 07 
3.445£: + 08 

1.321i; + 09 
5.069i; + 09 

1 945E + 10 
7 460f; + 10 

2 862f; + 11 

1 098 £: + 12 
4.212i; + 12 
1 61615 + 13 
6 198f; + 13 
2.40015 + 14 


0 00000 £ + 00 

6 18207£ - 10 
3 29844i; - 08 
1.96946£; - 07 
9 80469f; - 07 

3 54032£: - 06 
1.05575^-05 
2.98693£ - 05 

7 95623£' - 05 

1 97315i; - 04 

4.54910f; - 04 
9.64816£ - 04 
1.83881£ - 03 
3.311241: - 03 
5.601551: - 03 

8.778451: - 03 
1.387241: - 02 
2.101401: - 02 
3.078771: - 02 
4.185841: - 02 

4 969011: - 02 
5.530661: - 02 
6.528861: - 02 
8.559601: - 02 

2 012021 : - 01 


7 8611: + 00 
1 1391: + 01 

4 3231: + 01 
1 6401: + 02 
6 2201 : + 02 

2.3601: + 03 

8 9511: + 03 
3 3961: + 04 

1 2881: + 05 
4.8861: + 05 

1.8541: + 06 
7.0321: + 06 

2 667 E + 07 
1 0121: + 08 

3 8391: + 08 

1.4561: + 09 

5 5241: + 09 
2 0951: + 10 
7 9491: + 10 
3.0151: + 11 

1.1441: +12 
4.3391: +12 
1 6461: + 13 
6.2451: + 13 
2.4001: + 14 


0 oooooi: + 00 

1 285311: - 09 
1 507121: -08 
8.763101: - 08 
4 055111: - 07 

1 391811: -06 
4 029611: - 06 

1 12695£ - 05 

2 931901: - 05 
7.191161: -05 

1.617451: - 04 
3.369001: - 04 
6.38253£ - 04 
1.130221: - 03 
1.901951: - 03 

2 987551: - 03 
4 697581: - 03 
7.093281: - 03 
1.043221: - 02 

1 415971: - 02 

1.685821: - 02 
1.881061: - 02 

2 225861: - 02 
2 958931: - 02 
6.917271: - 02 


7.8631: + 00 
1.4831: + 01 
5.5531: + 01 
2.0801: + 02 

7 7S9E + 02 

2 9171: + 03 

1 0921: + 04 

4 0911: + 04 
1.5321: + 05 

5 737 E + 05 

2.149.E + 06 

8 0461: + 06 

3 0131: + 07 
1.1281: + 08 
4.2261: + 08 

1.5831: + 09 
5.9271: + 09 

2 2201: + 10 
8.3131: + 10 
3 1131 : + 11 

1 . 1661 : + 12 
4.3661: + 12 
1.6351: + 13 
6.1231: + 13 
2.4001: + 14 


sutron star mass of the radius and of the crustal thickness respectively are illustrated 
i Figs 4 and 5 From Fig. 4, it can be seen that for the EOS models Walecka (1974), 
id Sahu, Basu & Datta (1993), i.e. curves D and H, the slope for the lower values of the 
able mass is noticeably positive. This is because of the substantial contribution of 
tractive force brought on by the scalar <r-meson interaction present in these equations 
' state. It may be mentioned here that the c-interaction plays an important role at 
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intermediate densities providing nuclear saturation, and must be included in any field 
theoretical EOS model At high densities the contnbution of the vector meson (a;) 
interaction, which provides repulsive force dominates, and the slope becomes negative. 
The crustal thickness (A) is decided by how much mass is contained in the core and 
intenor regions of the neutron star, so that, in general, the more the mass of the neutron 
star, the less the value of A for any EOS 

As already emphasized, the crustal density profile of neutron stars is a key input in 
the proper understanding of a variety of important features associated with pulsars, 
such as the coohng rate Another important application, that is of much current interest 
in astrophysics, relates to the evolution of magnetic fields of pulsars However, several 
calculations reported in the literature on the problem of magnetic field evolution have 
considered cooling models for a specific EOS model while using the crustal model that 
corresponds to a different EOS model This is mainly because of non-availability of 
detailed, numencal crustal density profiles. The survey presented here is a preliminary 
step towards a uniform and consistent investigation of the problem of evolution of 
neutron star magnetic fields, which will be reported in a forthcoming paper (Bhat- 
tacharya & Datta 1995). 
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Abstract. Jupiter flux at 327 MHz was monitored using the Ooty radio 
telescope from July 12th to July 29th dunng the collision of comet 
Shoemaker-Levi 9 with Jupiter Flux was found to increase steadily from 
July 17th to July 26th by ~ 2-5 Jy, after which it declined to its pre-event 
value The comparison of 327 MHz observations with those at 840 MHz 
and 2240 MHz indicates that the enhancement was mainly due to the 
increased synchrotron emission and the contribution of thermal emission 
was very small at metric-decimetric frequencies The enhancement in radio 
emission was found to be more at 840 MHz than at 327 or 2240 MHz. The 
steepening of the spectrum between 327 and 840 MHz as well as between 
2240 and 840 MHz was also noted 

Key words SL9-Jupiter collision — radio flux monitoring — Jupiter radio 
emission 


1. Introduction 

Radio emission from Jupiter spans a wide frequency range from 001 MHz to 300 GHz 
This radio spectrum consists of three distinct components. Thermal radiation from the 
cloud layer dominates at millimeter and centimeter wavelengths At metric-decimetric 
wavelengths, the dominant component is the synchrotron emission by relativistic 
electrons trapped in the Jupiter’s magnetic field. At the low frequency end, below 
40 MHz, the emission consists of sporadic bursts due to plasma instabihties in the inner 
magnetosphere that generate the emission at frequenaes shghtly above the local 
electron cyclotron frequency 

In the metric-decimetric band of interest to us, the emission is due to synchrotron 
process powered by electrons accelerated to lO-lOOs of MeV by radial diffusion from 
the outer magnetosphere into ~ 2 Rj These high energy electrons are then trapped in 
a shell of ~ 2 J?j above the cloud top level and centered on the magnetic equator which 
is tilted by about 10° with respect to the spin equator. The muror points of most of these 
trapped electrons are within ~ 1 above and below the magnetic equatorial plane. 
The observed synchrotron source therefore is the projection of a toroidal source 
surrounding the planet. The width and height of this projection are approximately 6 Rj 
and 2 Rj respectively The spectrum in the metnc-decimetnc band is nearly flat with the 
oeak at around 1000 MHz. At higher frequencies, the thermal component dominates 
he synchrotron emission as shown in Fig. 1 
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wavelength 



Figure 1. Radio spectrum of J upiter observed over metric and decimetnc frequency bands The 
observed spectrum (solid line) is composed of emission due to thermal and non-thermal 
components (dashed lines), (Carr, T D etal 1986, Physics of Jovian Magnetosphere, pp 235). 


The passage of cometary fragments through the Jovian magnetosphere and the 

subsequent collisions were expected to affect Jupiter’s metnc-decimetnc emission. 

Vanous scenarios for the change m the Jovian radio emission were thought of and are 

mentioned below. 

■ It was expected that the energy released in the impact could ionize the cometary 
matenal resulting m an enhancement in the electron density These electrons could 
be accelerated to high energies in the shock waves that would have resulted from the' 
impact. The resulting increase m the number of relativistic electrons would enhance 
the emission. 

■ The emission could also be enhanced, if instead of producing new relativistic 
electrons in the impact, low energy electrons already existing m the synchrotron 
radiation belts were accelerated to sufficiently high energies thus resulting in a net 
enhancement in the density of high energy electrons. 

■ Enhanced thermal emission due to the heatmg up of the Jovian atmosphere by the 
impacts was also a possibility. Enhancement due to thermal emission was however 
expected to be seen for short durations after the impact because of short cooling time 
scales 

■ On the other hand, mcrease in the dust grams and larger-sized matenal in the 
magnetosphere with the entry of the comet was expected to slow down the high 
energy electrons As a result, the intensity was expected to decrease immediately after 
the cometary matenal entered the radiation belts Inward diffusion of electrons 
through a dusty magnetosphere was therefore expected to result in a drop m emission 
for many months after the comet had entered the magnetosphere (de Pater 1994). 
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We briefly descnbe our observational results in section 2 Results are discussed in 
section 3, and section 4 concludes the paper 


2. Observational results 

The observations were earned out using the Ooty Radio Telescope (ORT). The ORT 
consists of a 500 m long cylindncal parabolic antenna with an eSective collecting area 
of about 8000 sqm situated in southern India at a latitude of + 12°. The telescope 
operates at 327 MHz The telescope can be steered around its equatonal axis from 
— 4-5 to +4 5 hours which allows one to monitor a single source for a stretch of 
9 hours. In declination, the telescope can be pointed electronically from + 35 to — 35° 
by phasing the array of dipoles (Swarup et al 1971). The telescope can be operated 



Figure 2. Radio flux measurements carried out between 12th July and 29th July 1994 
Measurements at 327 MHz (a) were done using the Ooty radio telescope which we have 
compared with the reported observations at 840 MHz (b) and 2295 MHz (c) from Molonglo and 
J P L observatones, respectively All the flux values were normalized to a distance of 4 AU. 
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either in total power mode, where contribution from all the dipoles is added, or in 
phase-switched mode, where the signals from the north and south half of the telescope 
are correlated. In the phase-switched mode the ORT beam is 2 0° in the right ascension 
and 3' 5 sec{8) in declination 

The observations for the Jupiter-Comet collision were started on July 12th and 
continued till July 29th, 1994, well after the collision had taken place. Jupiter observa- 
tions were interspersed with observations of calibrators. A total of five calibration 
sources were used dunng the course of these observations, and at least three of them 
were observed on a given day (except for July 16th and 17th when the number was 2 and 
1 respectively). The data found contammated with interference were rejected. Observations 
were corrected for the effect of ORT beamshape as Jupiter dnfted in the beam For any 
Jupiter observation, confusion sources in the ORT beam down to 10 milli Jansky 
(1 Jansky = 10"^^ Wm“^Hz“^) were subtracted After calibrating it against the 
cahbration sources observed that day, a daily average of the Jupiter flux was obtained. 
The daily vanation of the Jupiter flux thus obtained is shown in Fig 2(a). 

As can be seen in Fig 2(a) there was a systematic increase of about 2-5 Jy (after 
normalizing it to a distance of 4 AU) observed between 17th and 26th July which was 
greater than the limit of flux measurement for ORT. This is indicated by the vertical 
error bar in Fig. 2(a) This increase in flux started after the impact of the first cometary 
fragment The maximum flux was observed on 26th July after whieh the flux steadily 
declined over the next 4 days of observations. The data from July 12th to July 16th, 
when the cometary fragments had entered the Jovian magnetosphere and were 
approaching Jupiter, show a gradual decrease in flux. 

We have compared our results with the observations at 840 MHz (Huntstead & 
Campbell-Wilson 1994) and 2240 MHz (Klein etal. 1994) (Fig 2b-c). The enhanced 
emission was more at 840 MHz than at the other two frequencies There is evidence of 
progressive steepening of the radio spectrum between 327 and 840 MHz, and between 
840 and 2240 MHz as a result of the impacts. As the radio emission declined after July 
22nd (at 840 MHz) the spectra in these two frequency ranges also became less steep. 


3. Discussion 

The initial decrease in the emission level (July 12th to July 16th) could have resulted 
from the absorption of high energy synchrotron electrons by the dust added to the 
atmosphere by the comet For the later enhanced radio emission we now explore two 
possible explanations -- both enhanced thermal as well as non-thermal emission. 

3 1 Thermal emission 

The data show a day-to-day increase in the flux with a peak on July 26th (i.e. 4 days 
after the last impact). The impacts therefore seemed to have had a cumulative effect on 
the radio emission rather than just a transitory enhancement at the time of the impact. 
The optical depth t for free-free absorption is (e.g. Benz 1993). 

T-0-llWe^v-^r''^(5L, (1) 

where is in number of electrons per cm^, v is in Hz, T in degree Kelvin and <5L is 
measured m cm. 
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In order to explain the enhanced emission at 327 MHz in terms of optically thick 
thermal emission, an enhanced electron density of the order of 10®cm"^ would be 
required, assuming T ^ 5000 K (since there is evidence that the temperature of the 
plumes resulting from the impact reached several thousand degrees, see Chapman 
1994) and 1000km. (Note that we cannot account for the observations by 
assuming optically thin thermal emission since the increase m flux at 840 MHz is more 
than that at 327 MHz) However gas with density as high as indicated above, will cool 
very fast. For example, taking an analogy with HI or HII regions since the dominant 
component of the Jovian atmosphere is hydrogen, the cooling rate can be roughly 
estimated by (Spitzer 1978) A/wJ ~7 x 10"^®ergcm^s“^ (assuming fraction of the 
H atoms ionized to be 10“^). 

The cooling time scale therefore would be U ^ {3/2nhkT)/A. Takmg nn 10^° (for 
fraction of ionized H atoms ~ 10‘^) and T 5000 K, we find that the cooling time 
would be ~ 1/2 hr Even if the cooling time derived above is off by an order of 
magnitude compared to what may hold for the Jovian atmosphere, we find that it is still 
difficult to understand the slow and sustained rise in the emission seen at the 
metnc-decimetnc band The therntal model however is more plausible for understand- 
ing the rapid variation in the emission reported at millimeter wavelengths 
Another difficulty with the thermal emission is that the electron density of 
10® cm" ® IS required across the Jupiter disk in a shell of thickness 1000 km. This is 

possible only if ALL the electrons from a typical cometary fragment spread out across 
the disk, which seems difficult to justify. 

We can however use the shape of the frequency spectrum to constrain the thermal 
contribution to the observed enhanced emission. Assuming the thermal contribution to 
be optically thick at 2240 MHz, and further assuming that the increase of flux of 
^ 1 56 Jy at 2240 MHz was purely due to thermal emission, in order to derive an upper 
limit on the thermal part of the enhanced emission, we find that the thermal contribu- 
tion at 327 MHz would be ^50mJy while at 840 MHz it would be smaller than 
300 mJy because the thermal emission scales as 1/2^ for optically thick emission 


3.2 Synchrotron emission 


The peak of the Jovian synchrotron spectrum is determined by Vmax(MHz) 4-8(E)^(B) 
sin(a), where energy and magnetic field is in MeV and Gauss respectively (Carr et al 
1986). Takmg v^ax 1 GHz, sin (a) 1 and 3 Gauss we find that the energies to 

which electrons need to be accelerated are 8 MeV to explain the enhanced emission. 
The power emitted by one electron is (Carr et al 1986): 


p«6’10 



sm^a^ 


3-10"'^W. 


( 2 ) 


Assuming that the enhanced emission comes from a volume of scale size 1000 km, 
we find that an excess electron density of the order of 2 5 x 10^ cm"® accelerated to 
energies 10 MeV is sufficient to explain the enhanced emission. 

Considering the moderate increase in the number density, compared to that for 
thermal emission, of relativistic electrons required to explain it, the synchrotron 
process seems to be the source of enhanced emission. The requirement of electrons 
accelerated to relativistic speeds could be met by a vanety of profcesses. 
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■ Shock waves generated by the impact could have accelerated the low energy 
electrons existing in the inner regions of the Jovian magnetosphere. 

■ There is an evidence of production of new electrons as inferred from the increase in 
auroras seen in ultraviolet images obtained by HST soon after the impact Therefore 
there is also a possibility that these newly generated electrons were accelerated and 
deposited in the inner regions of the magnetosphere. 

■ A third possibility would be the reconnection of magnetic field taking place behind 
the blob of ionized material as it arose after the impact and stretched the field hnes. 

All three of these scenarios would also explam the inward movement of the 
synchrotron belts. 


4. Conclusion 

We have observed the radio emission from Jupiter at the time of colhsion with comet 
Shoemaker-Levy at 327 MHz using the Ooty radio telescope. A careful analysis of the 
data indicates that there was a systematic but marginal increase of 2*5 Jy in the 
Jovian radio emission at 327 MHz due to the cometary impacts. The flux continued to 
mcrease till July 26th after which it declined. 

Comparison with observations earned out at 840 MHz and 2240 MHz indicates that 
the enhanced radio emission from Jupiter due to the cometary impacts was due to 
synchrotron emission. Acceleration of electrons in the mner regions of the magneto- 
sphere either due to shock waves emanating as an after effect of the impact or due to 
magnetic reconnection could be the cause of the mcrease in emission 
Though the thermal emission cannot explain the sustained rise m the radio emission 
at meter-decimeter wavelengths, it is more plausible for ‘bursts’ seen at the time of 
impacts as e.g in millimetric wavelengths. 
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Abstract The effect of intense magnetic field on the cosmic quark-hadron 
phase transition and also on the baryon number inhomogeneity has been 
investigated using phenomenological MIT bag model for the quark sector 
For the sake of simplicity an ideal gas equation of state has been considered 
for the hadronic phase 

Key words' Cosmic quark-hadron phase transition — Baryon number in- 
homogeneity — equation of state — Landau levels — quark nuggets. 


1. Introduction 

In the past few years a lot of work has been done on the effect of strong magnetic field on 
various properties of dense astrophysical and cosmological matter These works are 
mainly related to the emissions from magnetized neutron matter (Kammker etal 1991, 
f992), electron-induced electromagnetic and weak processes (ShuFman 1991), and the 
Landau diamagnetism and Pauli paramagnetism of dense nuclear matter (ShuFman 
1991, Fushiki et al. 1992) Very recently the effects of such intense magnetic field on the 
primordial nucleo-synthesis and also on the expansion rate of the universe micro- 
second after the Big Bang have been studied by Schramm et al (Cheng et all 993, 1 994) 
and Grasso & Rubinstein (1994) In recent years some interesting work has also been 
done on quantum electrodynamics in the presence of strong external magnetic field 
(Persson & Zeitlm 1994, Zeithn 1994; Danielsson & Grasso 1995) and also on the 
stability of matter in the presence of strong magnetic field (Vshivtsev & Serebryakova 
1994, Elliott et al 1995) In a senes of publications we have also reported some new 
results on the effect of strong magnetic field on first order quark-hadron phase 
transition at the core of a neutron star In particular, we have studied the stability of 
strange quark matter, equation of state of such stable phase and also the nucleation of 
quark bubbles in metastable neutron matter in the presence of strong magnetic field 
We have seen that strange quark matter becomes more stable m the presence of strong 
magnetic field and the equation of state of magnetized strange quark matter differs 
significantly from that of a non-magnetized one The most mteresting observation was 
that there cannot be any quark droplet nucleation in metastable neutron matter if the 
strength of the magnetic field is strong enough to make the quarks populate the Landau 
levels (Chakrabarty & Goyal 1994, Chakrabarty 1994, 1995a, b). 
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III all these reported results, mentioned above, the possible effects of strong magnetic 
fields on cosmic QCD phase transition have not been considered Initially the universe 
IS at a high temperature and in the quark-gluon plasma phase The net baryon number 
of the universe resides entirely in the pnmordial quark soup Due to cosmic expansion 
the temperature of the quark soup decreases and near the critical temperature, T^, 
a transition to colour singlet hadronic phase takes place During phase separation the 
baryon number is transported to the hadronic phase mainly by the formation of 
nucleons which are the lightest baryons. The other possible hadrons are mostly light 
mesons, which do not carry the baryon number Now, the quarks in the quark gluon 
plasma are almost massless and carry 1/3 baryon number, whereas the mass of the 
lightest baryon, the nucleon, is -- 940 MeV. A simple thermodynamic calculation 
shows that the baryon number solubility of pnmordial quark soup is much higher than 
that of hadronic matter at the critical temperature T^. This gives rise to inhomogeneity 
in the baryon number distribution This particular aspect has also not been studied in 
the presence of strong magnetic field. The aim of this paper is to investigate the effect of 
strong magnetic field on both cosmic quark-hadron phase transition and baryon 
mhomogeneity 

In § 2 we give the basic thermodynamic formalism of quark matter in the presence of 
strong magnetic field In sections 3 and 4 we study the effect of strong magnetic field on 
the cosmic QCD phase transition and on the baryon mhomogeneity Section 5 contains 
the conclusion of this work 


2. Formalism 

As is well known, the energy of a charged particle changes significantly in the quantum 
hmit if the magnetic field strength is equal to or greater than some critical value 

= mf l{q^h) in Gauss; where and are respectively the mass and charge 
(absolute value) of the particle (e g , = 2e/3 for u-quark, e/3 for d and s quarks, here 

e — \e\ is the absolute value of electronic charge), ft and c are respectively the reduced 
Planck constant and velocity of light, both of which along with the Boltzmann constant 
kg are taken to be umty in our choice of units For an electron of mass 0 5 MeV, the 
strength of this qntical field as mentioned above is B5J^^®^ = 44 x 10^ ^G, whereas 
for a light quark of current mass 5 MeV, this particular value becomes 10^ x 
arise for s-quark of current mass 150 MeV, it is ^ 10^® G, which is too high to realize 
at the core of a neutron star. But the possibility of such high magnetic field cannot 
be discarded in the case of the early universe (Vachaspati 1991). The critical magnetic 
field as defined above is the typical strength at which the cyclotron lines begin to occur, 
and m this limit the cyclotron quantum is of the order of or greater than the 
corresponding rest energy This is also equivalent to the requirement that the de- 
Broglie wavelength is of the order of or greater than the Larmor radius of the particle in 
the magnetic field. 

To study the cosmic QCD transition in the early universe m the presence of a strong 
magnetic field, we have considered the conventional MIT bag model. For the sake of 
simplicity we assume that quarks move freely within the system The current masses of 
both u and d-quarks are assumed to be extremely low (in our actual calculation we have 
taken the current mass for both of them to be 5 MeV, whereas for s-quark, the current 
mass is taken to be 150 MeV) 
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For a constant magnetic field along the z-axis (B„ = = constant), the 

single particle energy eigen value is given by 

+ mf + q,BJ2n + s + (1) 

where n = 0, 1, 2,. . being the principal quantum numbers for allowed Landau levels, 
s = ± 1 refers to spin up (+) or down (— ) states and k is the component of particle 
momentum along the direction of external magnetic field. Setting 2v = 2w + s + 1, 
where v = 0,1,2, . , we can rewrite the single particle energy eigen value in the 
following form 

= + ( 2 ) 

Now it IS very easy to show that the v = 0 state is singly degenerate while all other states 
with V 0 are doubly degenerate 

The general expression for the thermodynamic potential of the system at tempera- 
ture T = p~^ is given by 


-TlnZ = Y^,V + BV, (3) 

I 


where Z is the grant-partition function and the explicit form of the thermodynamic 
potential density for the ith species is given by 


n = -^ 

‘ ( 271 )" 


Jd^fcln(l -t-exp(J?(/i,-e,))), 


(4) 


where g^ is the degeneracy of the ith species ( = 6 for a quark or an antiquark), B is the 
bag pressure and V is the volume occupied by the system The sum in equation (3) is 
over M, d, s-quarks and their antiparticles We assume that for antiparticles the chemical 
potential /i, = — i.e they are in chemical equihbnum with respect to annihilation 
processes. 

Now let us consider the necessary changes to be made in equation (4) if a strong 
external magnetic field is present in the system Here we shall investigate the magnetism 
arising from quantization of orbital motion of charged particles in the presence of 
a strong magnetic field We know that if the magnetic field is along z-axis, the path of 
the charged particle will be a regular helix whose axis bes along the z-axis and whose 
projection on an x-y plane is a circle If the magnetic field is uniform, both the linear 
velocity along the field direction and the angular velocity m the x-y plane will be 
constant, the latter anses from the constant Lorentz force expenenced by the particle 
Quantum mechamcally the energy associated with the circular motion in the x-y plane 
is quantized in units of Iq^B^, The energy associated with the hnear motion along the 
z-axis IS also quantized; but in view of the smallness of the energy intervals, they may be 
taken as continuous variables. We thus have equation (1) or (2) as single particle energy 
eigen value. Now, these magnetized energy levels are degenerate because they result 
from an almost continuous set of zero field levels. All those levels for which the values of 
the quantity k\ + /cj he between and 2^jB„(v -1- 1) now coalesce into a single 

level characterized by the quantum number v. The number of these levels is given by 


S 

{Inf 



dk^dky = 


2n ’ 


(5) 
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here S is the area of the orbit m the x-y plane. This expression is independent of v. Then, 
in the integral of the form jd^kfik), we can replace jjdk^^dky by the expression given 
above, whereas the limit of k^, which is a continuous variable, ranges from - co to 
+ 00 . Then we can rewrite equation (4) for the thermodynamic potential density in the 
presence of a strong magnetic field in the form 

^ I I ” dkMi + exp(^(At, - et‘>))). (6) 

We have used equation (4) for 5 and s if < lO^^G and equation (6) for u, d-quarks and 
their anti-particles Whereas for equation (6) has been used for all the 

constituents In equation (4), = (k^ + mf ivhereas in equation (6), it is given by 

equation (2) 

The expression for kinetic pressure of the system is given by 

(7) 

I 

and the number density of the ith species is given by 



3. Cosmic quark-hadron phase transition 

Assuming for the sake of simplicity that only nucleons and pions are present in the 
hadronic phase and that they obey the free hadromc gas equation of state, the 
conditions to be satisfied by the co-existing phases are 

^ip = ^^d + 2ft,, (9) 

ftn = /*u + 2/id, (10) 

( 11 ) 

Ph = P,, ( 12 ) 

and 

(13) 

The chemical equilibrium conditions are also valid for anti-quarks and anti-hadrons 
i.e. we can write down equations similar to equations (9) and (10) for q, n and p (where 
q = M or d). We also assume that the condition of ^-equilibrium is satisfied in both the 
quark-gluon and hadronic phases For the anti-particles, the chemical potential 
—fii. The expressions for kinetic pressures are then obtained from equation (7) by 
replacing £1, either by equation (4) or (6), dependmg on whether or > 

respectively 

In Fig. 1 we have shown the phase diagram (n ^ vs T diagram) for such a mixed phase 
of quark-gluon plasma and hadronic matter. The dotted curve is for B„ = 0, whereas 
the solid curve is for = 10^ x This figure shows that the phase diagram 
changes significantly in the presence of a strong magnetic field. 
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Figure 1. Phase diagram (T—n^ plot) for a mixture of QGP and hadronic matter The dotted 
curve IS for = 0 and the solid one is for B = 10^ x 

in mm 


4 BaryoD inhomogeneity 

Before the quark-hadron phase transition in the early universe, the quark gluon plasma 
had almost an equal number of quarks and anti-quarks, since the baryon number of the 
universe is so small, 


(14) 

where n^ ad n^ are respectively the net baryon number and photon number densities in 
the early umverse This ratio is too small to affect the dynamics or the evolution of the 
early universe. Since the present universe is matter dommated, aU the anti-baryons 
must have been annihilated with their counterparts 
Now in the quark gluon phase, the baryon number is earned by the quarks, which 
are almost massless (except s-quark, whose mass is ~ 150 MeV) On the other hand, in 
the hadronic sector, the lightest baryons are nucleons of mass ~ 940 MeV. Therefore 
the baryon number prefers to reside m the quark gluon phase during phase transition 
(Witten 1984), or in other words, the baryon number solubility of primordial quark 
soup IS much higher than that of hadromc phase at We can estimate the ratio of 
baryon number densities in these two phases at T^, given by 


R = 



(15) 


% 

A simple non-relativistic thermodynamic calculation (Witten 1984) shows that 
R ~ exp(Mj/Tj), where is the hadronic mass » T^. Therefore the ratio J? is a large 
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Table 1. 



r^(MeV) 


00 

100 

455 4 


150 

33*4 


180 

14 7 


200 

99 

10" 

100 

2732 69 


150 

200 66 


180 

88 49 


200 

59 21 

10" 

100 

2783 01 


150 

201 39 


180 

8903 


200 

61 35 

10® 

100 

51908 2 


150 

1695 85 


180 

520 75 


200 

285 51 


number. It has been shown that « 6 or 100 for = 200 or 100 MeV respectively 
(Kurki-Suonio 1988, 1991; Fuller et al. 1988; Kapusta & Olive 1988) 

Since the net baryon number of the universe is small, fig = -► 0. Then the ratio R is 

given by 


R = lim 


1 i^^gld^>q)T = Tc 
3 


(16) 


Usmg the explicit form of 0^(1 = g or B), it is very easy to show that m the limiting 
condition n^-^O, the ratio R reduces to 0/0 form. Then by UHospital’s theorem, we 
have 


R = 


lim 


1 

3 


(17) 


which IS not in the 0/0 form. Using the explicit form of the second denvatives as written 
above we have calculated/? for four different values of T,(= 100, 150, 180and200MeV) 
and for = 0, 10^ x 10^ x and 10® x In Table 1 we have shown the 
variation of R with and . The variation of R with has already been reported and 
the results presented here are consistent with those of old published data. The vanation 
of R with has not been done before. If > 10^ x , R changes significantly from 
zero field values. If this is true for the early umverse, this could possibly lead to a very 
high baryon number inhomogeneity, which could modify the picture at the end of the 
quark-hadron phase transition. It is also obvious that the presence of a strong magnetic 
field favours the formation of quark nuggets dunng the early universe quark-hadron 
phase transition, which is expected to be one of the natural candidates for baryonic 
dark matter. 
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iVe have seen that in the presence of a strong magnetic field greater than the 
:or responding critical field, the phase diagram for a mixture of quark gluon plasma and 
ladronic matter differs significantly from that for zero magnetic field case. 

The other interesting observation is that the baryon number inhomogeneity in the 
larly universe increases by a few orders of magnitude if the magnetic field strength 
becomes greater than 10^ x 

If this theoretical prediction is correct, then it may have some effect on the primordial 
luclear abundances and also on the formation of quark nuggets during the first order 
luark hadron phase transition in the early umverse, which in turn affects the Big Bang 
lucleo-synthesis. 

Here we have given a very simplified picture of the baryon number transport during 
he first order cosmic quark hadron phase transition in the presence of a strong 
nagnetic field. Such calculations should only be used to indicate the order of 
Tiagnitude of the baryon number penetrability and a rough comparison with zero 
nagnetic field case. 
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Abstract. The UBV light curves obtained by Duerbeck (1975) and Hol 
(wide) and Ha (narrow) hght curves obtained by Chambliss & Davan (1987) 
of the detached echpsing binary VV Orionis (VV On) were analysed using 
the Wilson-Devinney method fixing the two parameters 7^(25, 000 K) and 
^(04172), resulting in the following absolute elements: A = 13.605 ± 
O.OSIRq, R* = 5.03 ± O.OSRq, R, = 2.43 ± O.OIR^, = - 5,18 ± 0.1 1, 
Afboi.c = — 1.54 ± 0.06, wifc = 10.81 ± 0.42mo and = 4.51 ± 0.41mp The 
de-reddened colours obtained from applying the reddening corrections of 
E{B~V) = 0”'.05 and £(C/-B) = 0'".04, and the denved temperatures of the 
components, gave spectral types of B1.5V for the primary and B4-5V with 
an UV excess of 0”. 3 for the secondary component. A companson of the 
log L and log of the components with the observed ZAMS shows the 
primary component to be a little above and the secondary component to be 
a little below/or on the ZAMS A comparison of the properties of the 
components of W Ori and a few other detached systems with the normal 
stars m the log L, log R and log versus log m planes, indicated a need for 
either a readjustment of the scales of the above parameters or modifications 
in the theoretical models. From the position of the components on the 
evolutionary tracks of Pop I composition computed by Schaller et al. (1992) it 
IS noticed that while the primary component of W On had shghtly evolved 
along the main-sequence, its secondary is still unevolved. The age of VV On 
IS found to be 10 ± 1 million years and it is at a distance of 368 ± 10 pc. 

Key words: Eclipsmg binaries — detached — W Ori — elements — mdividual. 


1. Introduction 

W Orioms (HR 1868 = HD 36695, P = i‘'.4854) was first found to be a light variable by 
Barr m 1903 (Barr 1905). It is an eclipsing binary of type EB with both the components 
belonging to the main sequence of spectral type B. There were many spectroscopic and 
photometnc investigations of this system. The spectroscopic studies of Struve 
& Luyten (1949) showed this system to be smgle-lined and they suggested a mass ratio 
of unity for the components. From their spectral studies, Beltrani & Galeotti (1970) 
estimated the mass ratio to be about 0.5. Using high dispersion spectra, Duerbeck 
(1975) could locate a few lines of the secondary component and was thus able to get an 
improved mass ratio of 0.45. However, the detection of the secondary lines had been 
questioned (Andersen 1976, Popper 1981) and the value of the mass-ratio remained 
ambiguous. 
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Fignre la. W On. Light curve in yellow. FiUed circles represent observations (normals) and 
the sohd line represents the theoretical curve obtained from the elements obtained from the 
combined V, B, Ha (n) solution (Table 2, col. 2). 



Figure lb. W Ori Same as m Fig. la, for blue 


Popper (1993) reported the detection of the lines of the secondary component of W On. 
His radial veloaty curves of the two components yielded a mass-ratio of 0.4171 Even 
though the variation of the center of mass vdodty Vq is poorly determmed. Popper (1993) 
suggests that “One can denve fundamental properties of W On, making use of the results 
m (his) Table 3, that should be improvements over previously published work”. Hence we 
felt it useful to re-analyse the published photometric data with the newly denved mass-ratio 
of 0.4172 and obtain reliable photometric elements for VV On The results of the 
analysis of a few published light curves of this system is given below. 


2. Earlier analysis of the light curves of VY Ori 

Of all the published light curves of W Ori, the UBV data of Duerbeck (1975) is the 
most extensive. Duerbeck himself had analysed his data usmg the Russell-Memll 
method (1952). Eaton (1975) analysed a series of OAO-2 light curves of this system 
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Figure Id. VV On Same as in Fig la, for ultraviolet Here, the theoretical curve is calculated 
using the elements given in (Table 2, col 3) 




AH<( W) 
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Figure le. W On. Same as in Fig la for Hot (w) Here, the theoretical curve is calculated usmg 
the elements given in (Table 2, col 4) 

which covered the wavelength range of 1300-4300 A, using a computational model 
developed by Lucy (1968) Cester et al (1978) analysed Duerbecks’s data using Wood’s 
WINK method. Budding & Najim (1980) analysed Eaton’s (1975) data with Kopal’s 
Frequency Domain method Chambliss & Leung (1982) solved their own UBVuvby 
hght curves with the Wilson-Devinney method. Giuricm et al. (1983) and Chambliss 
(1983) analysed the data of Chambliss & Leung (1982) with Wood’s WINK pro- 
gramme Chambliss & Davan (1987) analysed their Ha (wide) and Ha (narrow) band 
observations of this system using the WINK programme. All the above analyses, 
wherever needed, used a mass-ratio of 0.45. 


3. Present study 
3.1 Data 

Since an improved mass-ratio for W On (Popper 1993) is now available and as Wilson- 
Devinney (1971) synthetic light curve method (W -D) is regarded as the most reliable among 
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Table 1. W On. Elements obtained from the solution of individual hght curves using W-D 
method 


Element 

1 

Initial 

parameter 

2 

U 

3 

B 

4 

V 

5 

Ha (w) 

6 

Main) 

7 



25,000 

25,000 

25,000 

25,000 

25,000 

25,000 



15,550 

15,667 

15,544 

15,544 

15,560 

15,571 



04172 

04172 

04172 

04172 

04172 

0 4172 



860 

86.050 

86108 

86128 

86078 

86 108 

rh 

pole point 

0 368 

0.3700 

0 3630 

0 3640 

0 3572 

0.3612 


side back 


0.4100 

0 3989 

04006 

0 3901 

0 3962 




0.3846 

0 3765 

0 3777 

0 3698 

0 3744 




0.3971 

0 3879 

0 3893 

0 3804 

0 3855 

rc 

pole point 

0182 

01859 

01774 

01774 

01822 

01774 


side back 


01936 

01837 

01837 

01893 

01837 




01879 

01791 

01791 

01841 

01792 




01922 

01826 

01826 

01881 

01827 


U 

0.894 

0 9194 

0 9149 

09082 

0 8934 

0 9010 


B 

0 886 






L,IL,+ L, 

V 

0 878 







Ha(w) 

0901 







Ha(n) 

0903 







U 

0106 

0.0806 

00851 

00918 

01066 

0.0990 


B 

0114 







V 

0122 







H(x{w) 

0099 







Ha(n) 

0097 






/a 


00 

00 

00 

0.0 

00 

00 

*Xh 

U 

0320 







B 

0.390 







V 

0 325 







Ha (w) 

0.240 







Main) 

0.240 







U 

0 390 







B 

0.470 







V 

0 370 







Ha{w) 

0250 







Ha(n) 

0 240 






*A = ^c 


10 






*Gh = Gc 


0 25 







♦Fixed parameters 


all the available synthetic hght curve methods, we re-analysed the UBV hght curves of 
Duerbeck (1975) and Ha (w) and Hcc (n) light curves of Chambhss & Davan (1987) with 
this method usmg Popper’s mass-ratio of 0.4172 The data published by Duerbeck 
(1975) contains 655 17, 675B and 668 K observations of VV Ori As the memory of the 
computer available with the department is limited, to save computer time, we formed 
16317, 1675 and 171V normal points from these observations. During the phases of the 
primary eclipse, observations within a phase range of about 0.003 (1°) and at the 
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Table 2. VV On Elements obtained from the combined solution olB^V 8 l Ha (n) light curves 


Element 

1 

Combined 
solution for 

B, V&Ha(n) 

2 

Parameters used for fitting the U 
& Hoi (w) hght curves 

U Hol (tv) 

3 4 



25,000 

25,000 

25,000 



15,579 ±28 

15,579 

15,579 



04172 

0.4172 

04172 

r 


8609±005 

8609 

86.09 


pole point 

03582 + 0.0013 

0.3633 + 00016 

0 3554 + 0 0033 


side back 

0.3916 + 00012 

0 3994 + 00017 

0.3875 + 00034 



0 3709 + 00013 

0.3768 + 00016 

0 3677 + 0 0034 



0.3816 ±0.0012 

0 3883 + 00017 

0 3781+00034 


pole point 

01762 + 00012 

0 1838 + 00018 

01807 + 00034 


side back 

01823 + 00013 

0.1912 + 0 0018 

0.1875 + 00035 



01778 + 00012 

0 1858 + 00018 

0.1826 + 00036 



01813 + 00013 

0 1899 + 00019 

0.1864 + 00035 


B 

09077 + 00015 

0 9147 ±00017 

0 8938 ±00284 


V 

0 9023 + 0 0015 




Ha(n) 

09005 ±0 0270 




B 

0.0923 

0 0853 

0.1062 


V 

00977 




Ha{n) 

00995 



'3 


00 

0.0 

00 


*Pixed parameters 


remaining portions of the light curve points within a phase range of 001 (4°) were 
combined to form one normal point By using the normal points instead of the actual 
observations, the distortions in the hght curve, if any, might have been obliterated, but 
the final results will not be affected (Chambliss 1983) The Ha (w) and Ha (n) data of 
Chambliss & Davan (1987) contain 692 Ha (w) and 686 Ha (n) observations. All these 
observations were also formed into 209 Ha (w) and 213 Ha (n) normal points. For our 
analysis these normal points were used and the number of points in each normal was 
taken as its weight. These normal points are shown as filled circles in Fig. 1 (a-e). 

3 2 Methodology 

As already mentioned, the Wilson-Devmney (1971) method was used for solving the 
above light curves For initiating the W-D method, one needs reliable preliminary 
elements. For this purpose we used the values of the inclination, i, fractional radii, 
and and the luminosities, and for the UBV pass bands given by Duerbeck 
(1975). The luminosities for Ha (w) and Ha [n) pass bands were taken from Chambliss 
& Davan (1987), The values of temperatures and and the mass-ratio are from 
Popper (1993) The lunb-darkening coefficients for UBV pass bands were from Gryger 
(1972) and for Ha (w) and Ha (w) were from Chambhss & Davan (1987). These 
preliminary elements are given in Table 1, col. 2 
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Table 3a. W On Elements obtained from the UBV data of Duerbeck (1975) 


Element 

Method 


Duerbeck 

(1975) 

R-M 

Cester et al 
(1978) 
WINK 

Present studies 

W-D 




25,000 

25,000 

T„K 


— 

16,190 

15,580 



0450 ±0 014 

0 45 

0417 

f 


861±02 

863±05 

8609 + 005 



0,368 + 0002 

0364 + 0003 

0 3702 ±0 0012 



0182± 0 002 

0.177 + 0002 

01784 ±0.0012 

h 

u 

0 894 

0.952 

0915 


B 

0 886 

0.968 

0.908 


V 

0.878 

0.965 

0902 

L, 

u 

0106 

0.048 

0085 


B 

0114 

0.032 

0092 


V 

0122 

0035 

0098 

mJniQ 


760 

7.6 

10 81 

mJniQ 


3 42 

3.4 

4 51 

RJRo 


447 

4.44 

5 03 

Rc/Ro 


2 21 

2.13 

243 

A/*,. 


-2 58 

— 

-2 59 

r{pc) 


372 

— 

368 


♦Assumed parameter. 


Table 3b. W On Elements obtained from the Ha (w) & Ha (n) data of Chambliss & Davan (1987). 


Element 

Method 

Wavelength 

Chambliss & Davan (1987) 
WINK 

Present studies (individual solutions) 
W-D 

Ha(w) 

Ha(n) 

Ha(w) 

Ha (n) 


25,000 

25,000 

25,000 

25,000 

Te,cK 

15,690 ±150 

15,760 ±170 

15,560 ±50 

15,571 ±54 


045 

0 45 

04172 

04172 

f 

83 6 + 0.5 

85.5 ±0 8 

86 078 ±0140 

86.108 ±0150 

h 

0.367 + 0.003 

0.361 ±0.004 

0.3691 ±0.0035 

03737 ±0 0036 

rc 

0179 ±0.003 

0.173 + 0.004 

01848 ±00036 

0.1798 ±00038 

L, 

0901 

0903 

0893 

0901 

Lc 

0 099 

0097 

0107 

0099 

h 

00 

00 

00 

00 


♦Assumed parameter 


According to the principles of the W-D method we adjusted the following par- 
ameters. the inclination i, the surface potentials £i^ and the relative monochromatic 
luminosity L^,, the temperature of the secondary and the third light The 

temperature of the pnmary component, 7^, ^^ (25,000 K) and the mass-ratio ^(04172) 
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Table 3c. VV On Elements obtained from the UBVuhvy data of Chambliss 
& Leung (1982) 



Chambliss & 

Giuncin et al. 

Chambliss 

Element 

Leung (1982) 

(1983) 

(1983) 

Method 

W-D 

WINK 

WINK 

*Te.HK 

21,000 

26,000 

25,000 

T,,K 

12,835 

17,017 

15,957 


04568 + 0 0012 

040 

045 


90 + 0.49 

86 3 

85 60±017 

h 

03530 + 0.0004 

0362±0003 

03627 + 00017 


0 1709 + 0 0004 

0.177 ±0004 

01764 + 00014 

L, 

U 09446 

0924 

0 925 


B 09115 

0907 

0 895 


V 09026 

0904 

0 891 

k 

U 00554 

0076 

0075 


B 00885 

0093 

0105 


V 00974 

0096 

0109 


0.0 

00 

00 

*Assumed parameter 



Table 3d. 

VV On. Elements obtained from the OAO-2 data 


Element 

Eaton (1975) 

Buddmg & Najim (1980) 

Method 

Computational method of Lucy (1968) Kopal’s frequency domain 


21,000 


— 

V 

1 1,800 ±500 


— 


0458 ±0050 


— 

p 

845 + 05 


87 5 ± 1.3 


03790 + 0.0040 


0.3708 + 0.0026 


01745±00040 


01748 + 00080 


* Assumed parameters 


were treated as fixed parameters. It is expected that analysis of the hght curves with 
fixed Ti and q at known values should yield most reliable elements The limb-darkening 
coeffiaents and the albedos and A^. and the gravity darkening coeflBcients 
and Gg were also kept constant. From our earlier studies we noticed that the choice of 
these coefficients had negligible effect on the derived elements. As it is now well 
established that VV Ori is a detached system, we used code-2 of the W-D method 
meant for detached systems. A sufficient number of runs of the DC programme was 
made till the sum of the residuals (0-C)^ showed a minimum and the corrections to 

the parameters became smaller than their probable errors In order to check the 
internal consistency of the results (Popper 1984) separate solutions for each of the 
UBVHa (w) and Hol («) light curves were made. These results are given m cols. 3, 4, 5, 
6 and 7 of Table 1. In none of the solutions was the third light present. 
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From the results of the above analysis, it is noticed that the relative fractional sizes and 
denved from the U solution, are slightly larger than those derived from B, V and Ha. (n) 
pass bands. Since only a small discrepancy exists between the and r, values denved from 
Lhe Ha (w) and Ha (n) pass bands, and that the and values of Ha (n) solution are almost 
identical with those of B and V solutions and m order not to combine many light curves 
and comphcate the analysis, it was decided to get a combined solution only for B, V and Ha 
;n) and as before, keeping T^ ,, and q as fixed parameters, another solution was made for the 
combined B, V and Ha (n) light curves. These results are given in Table 2, col. 2 The third 
li gh t /j IS absent in this solution also. Using these elements, we have fitted the U and Ha (w) 
light curves. In doing so we had kept all the parameters except as fixed parameters The 
final parameters used for fitting the U and Ha (w) light curves are given in Table 2, 
cols 3 and 4 For calculating the absolute elements, we used the elements obtained from 
the combined solution only However for obtaining the {U-B) colours of the compo- 
nents, we used the ^ and ,, values given in Table 2, col 3 
The theoretical curves obtained from the elements given in Table 2 for the UBVHa (w) 
and Ha (n) light curves are shown as solid lines in Fig 1 (a-e) The fit of the theoretical 
curves to the observed normals is quite satisfactory. In Table 3a we have given for 
comparison the elements denved from the UBV data of Duerbeck (1975) Table 3b gives 
the elements obtained by Chambliss & Davan (1987) from the analysis of their Ha (w) and 
Ha (n) data In Tables 3c and 3d we have compiled the elements obtained by various 
mvestigators using the same data but different methods of solution It can be noticed that 
the values of the denved parameters depend only on the method of analysis and not on the 
underlying data. In such analysis the most significant differences exist between the 
parameter i and L The high value of 90° obtained for i from the UBVubvy data usmg the 
W-D method by Chambliss & Leung (1982) (Table 3c) may be attributed to their 
combining many light curves which might have complicated the solution. 


4. Spectral types 

From the denved lummosities of V and B pass bands for the combined solution (Table 2, 
col. 2) and U pass band (Table 2, col 3) and usmg the corresponding differential 
magnitudes for unit luminosity at the quadrature [AK= — 0"'885, AB= — 0'"902, 
AU = — ('".OOfi] and the standard magnitudes and colours of F = 6™.24, {B-V) = 
- 0”’.18 and (U-B) = - 0'".8 1 for the comparison, Cj , (Duerbeck 1975) we obtained the 
following magnitudes and colours for the individual components of VV On; 



Hot 

Cool 


component 

component 

V 

5.46 

7.88 

(B-V) 

-021 

-014 

iU-B) 

-092 

-0.82 


According to Eaton (1975) an interstellar reddening correction of E(B-V) = O'" 05 is 
apphcable to VV On. Assuming the law of interstellar reddening (Allen 1976) to be 
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holding good, this reddening in E(B-K) corresponds to a reddening of E(U-B) = 0 04 
and /4„ = O'”. 16, Correcting the above derived magnitudes and colours for the redden- 
ing, one gets 



Hot 

Cool 


component 

component 

Vo 

5.30 

7 72 

(B-V), 

-0.26 

-019 

iU-B)o 

-0.96 

-0.86 


The above corrected (B-K)o of primary component corresponds to a spectral type of 
either B2 V (Allen 1976) or B1 V (Popper 1980), giving it an average spectral type 
B1.5 y. For such a star the (U-B)o would be — 0 96 (Allen 1976; Schmidt'Kaler 1982) 
which IS the same as that obtained above. Hence both the derived (B-POo (^"B)o 
colours give a consistent spectral type of B 1.5 F to the primary component The 
average temperature of a B1 5 F star would be 25,500 ± 800 K (Allen 1976, Popper 
1980; Schmidt-Kaler 1982), which is the same as the fixed parameter, ^ (25,000 K) 
used for solvmg the light curves. This confirms that our analysis with fixed 7), (25,000 K) 
and ^(0 4172) yielded results akin to the reddemng free data. 

The denved (B-y)o of - 0".19 of the secondary component corresponds to a spectral 
type of B4 V (Men 1976), and its denved temperature of 15,580 K suggests a spectral 
type of B5 V (Allen 1976, Popper 1980). Thus both the (B-POo colour and temperature 
give a consistent spectral type of B4-5P^ to the secondary component. However, for 
such a star, the (U-B)o should be -0'”.60 (Allen 1976, Schimdt-Kaler 1982), but its 
denved (l/-fl)o = — 0”. 86 Hence the secondary component appears to be a B4-5P^ star 
with a LTV' excess of about 0” 3 In view of their small sizes and low luminosities (section 
5) the possibility of the colours of these components corresponding to those of super- 
giants is not considered According to Popper (1993) the properties of the components 
of W On correspond well with those of the main-sequence stars of types B1 and B5. 
Eaton (1975) classified the primary as BO V to B1.5 V while the secondary in the middle 
B range and concludes that ‘VV On seems to be a reasonably normal mam-sequence 
pair’. Chambliss (1983) suggested that W OnA is a B1 K star Duerbeck (1975) derived 
a spectral type of B1 V for the primary and B5-9 V for the secondary components. 
Hence our derived spectral types of the components (Bl,5 V + B4-5 V) agree quite well 
with those of other investigators except that the secondary appears to have aUV excess 
of 0”.3. Since the standard magnitudes and colours of the comparison star used by 
Chambliss & Davan (1987) are not available, we could not get the magnitudes and 
colours of the components of VV On in Hct pass bands. 


5. Absolute elements 

Taking the amplitudes of the radial velocity curves as = 136 ± 6 and K 2 = 326 ± 
7 km/sec from Popper (1993) and i = 86.09^ from the present study (Table 2, col. 2) one 
gets the masses of the individual components of VV On as = 10.81 ±0.42 and 
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mJmQ = 4.51 ± 0 41 and A/R^ = 13.605 ± 0 031 Combining these values with those of 
r^, Tf^ and we obtained the following sizes, luminosities, bolometnc and visual 
magnitudes 



Hot 

component 

Cool 

component 

R/Ro 

5.03 ±003 

2.43 ±0.02 

logL/Lo 

3.948 ± 0044 

2.493 ±0023 


-5.18 ±0.11 

-1.54 ±0.06 


-2.59 ±0.11 

-0.10 ±0.06 


In calculating we used the bolometric correction of — 2"* 59 and — l”44forthehot 

and cool components respectively (Popper 1980). Applying the interstellar absorption 
of = O'". 16, the distance of W On was found to be 368 + 10 pc 


6. Component properties and comparison with other similar systems 

From a plot of the positions of the two components on observed ZAMS in Fig 2 (Morton 
& Adams 1968) it is noticed that the primary component hes slightly above and the 
secondary component slightly below or on the ZAMS. Chambliss (1984) found the larger 
component (primary) to be normal for a spectral type of Bl V and found it to be fairly close 
to the ZAMS Eaton (1975) concludes that W On seems to be a reasonably normal 
mam-sequence pair and when plotted on the observational ZAMS the primary lies 
a little above the main-sequence and the secondary lies slightly below. 

In order to compare the above properties of VV On with other detached systems 
havmg B type components, we plotted in Fig 2 the positions of a few systems as given in 
tables 2 and 4 of Popper (1980) From this we notice that except a few primanes, 
notably no. 14 (%^ Hya), all the other components he either on the ZAMS or very near to 
it. A little adjustment m the scale of either the denved parameters or the observed 
ZAMS in the log range of 4.3 to 4 05 may make most of the systems to fall on the 
ZAMS. For ready reference, the list of systems given by Popper (1980) and used in this 
study are given in Table 4 

In Figs 3, 4 and 5 (solid hnes) we have plotted the relations of log L, log R and log Tg 
with respect to log mass of the normal mam-sequence stars (Allen 1976) In these 
figures, in addition to the positions of the components of VV On, we have also plotted 
the positions of other detached systems given in Table 4. From these plots one can see 
that the primary component of VV On is overluminous, hotter and a little undersized 
for its mass. In the case of the secondary component it is seen that while its luminosity is 
consistent with its mass, its radius is smaller and the temperature is higher for its mass. 
However from a companson of the positions of the other detached systems m these 
plots, It IS apparent that in addition to gettmg more accurate elements, there is a need 
for either refixmg the scales of the observed log L, logR and log Tg or modifying the 
models and gettmg new theoretical relations. Unless this is done, one cannot come to 
any conclusion regarding the properties of the components of VV On and compare 
them with other detached systems. 
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4.5 4-4 4.3 4-2 4-1 4-0 3-9 


log Te 

Figure 2. W On Plot of log L/Lq versus log Tg. The solid line represents the observed ZAMS 
(Morton & Adams 1968). The primary (P) and secondary (5) components of VV On are indicated 
as ‘-I-’ Components of other detached systems are shown as ‘x* (P) and ‘A’ (5) Numbers 
represent the names of the systems given in Table 4 ‘®’ represent both the pnmary and 
secondary components of a system 


7. Evolution 

With a view to study the evolutionary status of the components of VV On, we used the 
evolutionary tracks computed by Schaller et al (1992) for stars of Pop I composition. 
These evolutionary tracks for stars of initial masses of 10 8 and 4.5 m solar umts are 
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ible 4. VV On List of detached systems used m Fig. 2 

3 Name logwii logL^ logi?i log logiWj logL 2 log /?2 log ^2 


YCyg 
V478 Cyg 
V453 Cyg 
CWCep 
I a Vir 
. QX Car 
' V539 Ara 

! CV Vel 
> DI Her 
) UOph 
L V760 SCO 
I AG Per 
5 ^Phe 
X Hya 
5 AS Cam 
5 V451 Oph 

7 RXHer 
3 AR Aur 


1.220 

445 

1 190 

4 62 

1.160 

4 69 

1070 

428 

1030 

4 33 

0960 

3 59 

0 790 

3 26 

0 780 

319 

0 720 

2 65 

0 710 

290 

0 700 

271 

0.660 

2 89 

0 593 

2 50 

0 558 

249 

0 520 

217 

0444 

186 

0 439 

1.79 

0.394 

167 


0.780 

4.485 

0.860 

4.485 

0.930 

4 470 

0 730 

4 470 

0910 

4 390 

0 630 

4.340 

0 640 

4 255 

0 610 

4 255 

0410 

4.220 

0 540 

4.220 

0.480 

4 200 

0.480 

4 245 

0455 

4160 

0 642 

4.063 

0 431 

4088 

0.423 

4.015 

0 387 

4.015 

0.262 

4048 


1220 

4.45 

1 190 

4 62 

1050 

4 20 

1040 

413 

0 830 

316 

0 930 

3.48 

0 720 

3 02 

0.780 

319 

0.660 

2.52 

0.660 

2 74 

0.670 

2 53 

0 610 

2 62 

0 405 

151 

0 422 

164 

0 400 

156 

0 373 

154 

0 367 

148 

0 360 

146 


0 780 

4485 

0 860 

4 485 

0 730 

4 445 

0.700 

4 445 

0 640 

4 230 

0 600 

4 330 

0.570 

4 230 

0.610 

4 255 

0410 

4190 

0.490 

4 200 

0 430 

418 

0410 

4 210 

0 267 

4005 

0 334 

4.005 

0 301 

4.001 

0 326 

3.985 

0 292 

3 985 

0.262 

3 997 


liown in Fig. 6, which is a plot of log LjL^ versus T^. In the same diagram, the 
iochrone for an age of 10 million years is also shown This isochrone was obtained by 
iterpolation between the evolutionary tracks of Schaller et al. (1992) based on points 
if corresponding evolutionary status One can see from this figure that the primary has 
volved from its initial position along the IO-ShIq track The secondary component is 
till unevolved and lies at the foot of the track of a 4 5 iMq star with a slightly lower 
uminosity. Both the components satisfactorily fit the isochrone of lOmilhon years and 
lence we conclude that the age of VV Ori is 10 i 1 milhon years 


8. Discussion 

rhe VV excess of about found for the secondary component may be explained 
IS due to hot gases surrounding this component, the source of which might be the gases 
ixpelled by the primary during its present evolution. Comparing the observed frac- 
lonal radii of 0.370 and 0.178 of the primary and secondary components with their 
intical Roche lobes of 0.457 and 0.299 respectively, (for a mass-ratio of 0.42) (Plavec 
fc Kratochvil 1964), it is found that these components are filling 80% and 60% of then- 
respective Roche lobes. Even though any stellar activity hke gas streams etc., are not 
generally expected in detached systems, a sUghtly bigger fractional radius (of about 
D.Ol) obtained for U pass band in our analysis suggests some form of gases or gaseous 
streams in the system From their studies of W On light curves, Duerbeck (1975) and 
Giuncin et al. (1983) had also arnved at similar conclusions. 

Even though the geometrical and absolute elements of VV Ori are now well 
established, a few peculiarities are still left unexplained. The infiuence of mutual 
interactions on the evolution of the system should be carefully estimated. 




IS log m/mQ The solid line represents the relation for nonnal stars (Allen 1976) The pnmary (P) and 
shown in the figure with error bars Other details are the same as for Fig 2 
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Figure 5. W On Plot of log versus \ogm/m^ Details are the same as for Fig 3 
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4-5 4-4 4-3 4-2 4-1 4-0 3-9 3-8 3-7 3*6 

log T<2 


Figure 6. VV On Evolutionary tracks of stars of initial masses 10.8 and 4 5 in solar units. These 
tracks were computed by Schaller et al (1992) for stars of Pop I composition The isochrone for 
10 X lO^yrs is also shown in the diagram, the positions of the pnmary (P) and secondary (S) 
components of VV On are shown as ‘ ’ 
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